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SECTION  1 
IHTRODOCTIOH 


The  progress  an£  result*  Achieved  on  this  basic  research  program  are 
sumsrized  herein.  The  purpose  of  the  research  program  is  to  study  the 
applicatiot.  of  near-field  techniques  to  charccteriae  the  radiation  and 
coupling  characteriacica  cf  wideband  CW  or  pulsed  radirtors.  Out-of-band 
as  well  as  in-band  situations  are  of  interest  and  both  situations  are 
included  in  the  study. 

Tne  radiation  pattern  of  an  antenna  can,  in  principle,  be  calculated 
by  conventional  deterministic  analysis  based  on  a  knowledge  of  all  system 
variables.  The  feed  netirark,  however,  can  support  multi-mode  energy 
propagation  at  out-of-band  frequencies.  For  example,  energy  can  be 
propagated  in  the  18  different  modes  delineated  in  Figure  1-1  in 
standard  Wk-284  S-band  waveguide  at  the  out-of-band  frequency  of  10>0 
GHz.  Furthermore,  all  of  the  modes  whose  second  index  is  non-zero  have 

an  electric  field  component  polarized  orthogonal  to  the  dominant  in- 

band  TE^q  mode.  The  electric  fields  for  the  first  5  modes  are  depicted 
schematically  in  Figure  1-2. 

Celculation  of  the  relative  phases  end  amplitudes  of  the  different 
modes  requires  a  detailed  knowledge  of  the  feed  network.  This  is 
generally  a  very  difficult  boundary  value  problem.  Further,  the 

coefficients  of  the  modes  will  be  quite  sensitive  to  minor  electrical 
and  mechanical  variations  in  the  feed  network,  especially  those 
containing  active  devices,  sc  that  nominally  identical  systems  can 

produce  distinctly  different  out-of-band  patterns.  This  characteristic 
out-of-band  pattern  sensitivity  is  evident  from  experimental  and 
theoretical  studies  [  1-7]  .  In  order  to  account  for  these  seemingly 
random  effects,  the  cut-of-band  mode  excitation  coefficients  may  be 
treated  as  random  variables.  Treatsient  of  these  stodal  coefficients  as 
randoo  variables  serves  a  dual  purpose.  First,  it  allows  the  analysis 
of  radiation  patterns  to  proceed  iodependent  of  the  feed  structure 
(which  may  vary  between  otherwise  identical  anteonas).  Second,  it 
lead*  to  statistical  equations  that  can  account  for  variations  in 
fabrication  and  excitation  of  the  nominal  antenna. 
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Figure  1-1.  Uagram  depicting  the  allowed  higher  order  aodes  versus 
frequency  for  WR- 284  S-bond  rectangular  waveguide. 


(c)  TEoi  MODE 


td)  TE„  MODE 


(e)TM„  MODE 


Figure  1-2.  Sketches  depicting  the  transverse  electric  field’  for  the 
indicated  higher-order  nodes  in  rectangular  waveguide. 
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Consequently,  ststistical  analysis  techniques  are  needed  for  out- 
of-band  characterization  since  it  is  iapractical,  or  perhaps 
impossible,  to  solve  the  electromagnetic  boundary  value  problems  for 
the  aorriad  reflector  and  phased  array  feed  systems  encountered  in 
practice.  Statistiral  techniques  and  concepts  can  also  be  employed  to 
provide  succinct  EMC  descriptions  of  essentially  "deterministic"  in- 
band  wideband  or  pulsed  radiators.  Consequently,  a  considerable  effort 
hac  been  devoted  to  statistical  characterization  of  wideband  radiators. 

The  research  progran  is  divided  into  four  major  tasks.  The  four 
tasks  are  as  follows. 

Tisk  1.  Provide  a  near-field  methodology  to  characterise 
electrosugnetic  emitter  radiation  patterns  at  in-band  and  out-of-band 
frequencies  for  wide  bandwidth  radiators.  The  objective  cf  this  task 
is  to  develop  the  appropriate  theory  and  equations  based  on  statistical 
analysis  techniques  for  efficient  characterization  of  wideband 
radiators. 

The  objective  of  Task  I  was  achieved  via  the  following 
accomplishments: 

(a)  The  theory  and  equ.'.tions  were  developed  for  both 
in-band  and  out-of-band  frequencies  for  deriving 
statistical  average  far-ficld  patterns  from  wide-band 
or  pulsed  near-field  measurements. 

(b)  Based  upon  the  theoretical  studies,  methods  of  efficiently 
characterizing  both  wideband  continuous-wave  and  pulsed 
radiators  using  near-fisld  measurement  Cecbniquea  were 
devised. 

(c)  Based  upon  the  theoretical  studies,  a  mnserical  jisKilation 
whicn  is  applicable  to  various  antenna  types,  such  as 
reflector  and  phssed-array  antenn4.s  was  performed  to 
demonstrate  that  valid  statistical  far-field  pattern 
distributions  can  be  obtained  from  near-field  measurements. 

Task  2.  Theoretically  relate  the  radiation  pattern 

characterization  to  the  basic  data  needed  for  efficient  optimization 
analysis  of  electromagnetic  spectrum  ucage.  The  objective  of  this  task 
is  to  relate  the  near-field  derived  wideband  antenna  characterisation 
to  antenna  coupling  of  antenna  systems  which  co-exist  in  the  sasm  EM 
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envivoraBeac .  This  objective  wes  achieved  by  the  following 
accoaplishaents : 

(a)  The  theory  and  equations  were  developed  for  analysing 
the  coupling  between  co-sited  transBitting  and 
receiving  ayateaw  for  co^inationa  of  wideband 
transnitting  antennas,  narrowband-haraonic  transaitting 
antennas,  and  both  wideband  and  narrowbacd'-bamonic 
receiving  antennas. 

(b)  Baaed  upon  the  theoretical  studies,  the  near-field 
derived  data  required  for  efficient  EM  spectrum  usage 
optimization  were  defined  and  employed  in 
numerical  simulations  to  demonstrate  application  to 
EM  spectrum  usage  analysis. 

Task  3.  Provide  the  methodology  to  assess  the  effects  of  system 
devices  (i.e.,  higher-order  mode  generation)  on  the  radiation  pattern. 

The  objectives  of  this  task  are  (])  to  determine  e  method  whereby 
the  pattern  effects  of  higher-order  modes  which  are  generated  by  system 
devices  at  out-of-band  frequencies  can  be  assessed  and  (2)  to  study  the 
feasibility  of  identifying  the  aK>dal  content  of  at.  out-of-bard  feed 
system  from  measured  d^ta.  These  objectives  were  achieved  by  the 
following  accomplishments. 

(a)  Equations  were  derived  during  Task  1  and  Task  2 
to  permit  the  out-of-band  radiation  pattern 
statistics  to  be  computed  from  a  knowledge  of 
the  system  device  statistics  which  describe  the 
higher-order  mode  generation  and  propagation. 

(b)  Theory  and  equations  for  identifying  the  higher-order 
mode  statistics  of  system  devices  through 
utilization  of  wideband  or  pulsed  near-field 
measurements  were  derived. 

Because  the  first  objective  for  this  task  was  essentially  achieved 
me  a  result  of  the  work  on  Tasks  1  and  2,  the  scope  of  the  effort  on 
this  task  was  ezpanded  to  include  a  brief  analytical  study  of  the  out- 
of-band  characteristics  cf  three  cosawn  waveguide  components.  The 
three  waveguide  devices  studied  are  (1)  a  coaa-to-waveguide  adapter, 
(2)  radial  bends,  and  (3)  a  ferrite  phase  shifter. 


T«8k  4.  Investigate  the  i^act  of  site  effects  on  the  near-field 
antenna  analysis  technology. 

The  objective  of  this  task  is  to  extend  the  existing  uunochrovatic 
spectruB  scattering  aatrix  analysis  to  study  antenna  siting  effects  on 
the  wideband  and  out-of-band  perforaance  of  radiating  systens.  This 
objective  was  achieved  by  the  following  acconplishnents. 

(a)  A  theoretical  study  was  pcrfonaed  to  investigate 
techniques  and  derive  the  equations  for  extending 
the  existing  Bonochroaatic  spectrum  scattering 
mat'-ix  theory  to  analyxe  site  effects  over 

wide  frequency  bandwidths. 

(b)  Based  upon  the  results  of  the  theoretical  study, 
a  numerical  simulation  of  the  site  effects  on  the 
wideband  in-band  and  out-of-band  performance 

of  a  radiating  system  was  performed. 

All  of  the  objectives  set  forth  in  the  task  statements  were 
successfully  achieved.  The  salient  results  of  the  research  work  for 
each  task  are  successively  presented  and  discussed  in  Sections  II 
through  V.  Host  of  the  important  results  for  Task  1  were  previously 
presented  in  the  Interim  Technical  Reports  Mo.  1  and  Ho.  2  [  4,5]. 
Similarly,  the  main  results  for  Task  2  «'ere  included  in  Interim 
Technical  Report  Mo.  2.  Abbreviated  versions  of  the  previously 
reported  results  for  Tasks  1  and  2  are  contained  in  Section  XI  and 
Section  III,  respectively,  along  with  discusaions  of  new  or  additional 
results.  Results  for  Tasks  3  and  4  are  contained  in  Sections  IV  and  V, 
respectively.  Conclusions  and  recoBMndations  baaed  on  the  entire 
research  progiam  are  presented  in  Section  VI. 
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SECTION  II 

THEORY  AND  TECHNIQUES  FOR  WIDEBAND 
ANTENNA  PATTERN  CHARACTERI2ATIGN 


A.  Introduction 

The  theory  and  equations  were  developed  for  characterizing  the 
radiation  patterns  of  wideband  cw  or  pulsed  antennas  over  both  in-band 
and  out-of-band  frequency  intervals  from  measured  data  collected  via 
near-field  measurement  techniques.  The  results  are  applicable  to 
either  phased  array  or  reflector  antennas.  Numerical  simulations  were 
performed  for  (1)  a  20-element  out-of-band  waveguide  phased  array  with 
no  inter-element  coupling  and  (2)  3-element  and  9-element  wire  arrays 
with  inter-element  coupling.  The  effect  of  statistical  correlations  of 
the  near-field  data  was  studied  and  methods  for  handling  correlation 
effects  were  derived.  Useful  approximations  for  the  probability 
density  function  for  the  radiated  power  pattern  statistics  for 
correlated  sources  were  also  identified. 

The  key  results  obtained  via  the  Task  1  efforts  may  be  summarized 
as  follows: 

(1)  The  statistical  average  patterns  and  standard  deviations  at 
selected  frequencies  can  provide  a  very  succinct  engineering 
description  of  the  important  EMC  characteristics  of  wideband 
CW  multimoding  antennas.  The  statistical  average  patterns 
and  standard  deviations  are  a  practical  alternative  to  the 
comparatively  more  expensive  and  cumbersome  Monte  Carlo 

s Imulations . 

(2)  The  statistical  average  pattern  for  a  given  frequency  may 
be  computed  from  a  knowledge  of  the  following  near-field 
statistical  parameters: 

(a)  statistical  average  value  of  the  electric  field  at  all 
near-field  measurement  points, 

(b)  the  standard  deviation  of  the  electric  field  at  all 
measurement  points,  and 

(c)  the  covariance  functions  for  the  electric  fields  at 
all  different  near-field  measurement  points. 

(3)  The  statistical  average  pattern  versus  time  for  a  pulsed 
system  depends  on  all  of  the  above  near-field  statistical 
parameters  listed  in  Item  2,  and  the  following  far-field 
statistical  parameters: 
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(a)  Che  sCaciscical  average  value  of  Che  far-field 
elecCric  field  at  all  frequencies  in  the  frequency 
band , 

(b)  the  standard  deviation  of  the  far-field  electric 
field  at  all  frequencies  in  the  frequency  band,  and 

(c)  the  covariance  functions  of  the  electric  field  at 
all  different  frequencies  in  the  frequency  band. 

However,  the  far-field  statistical  parameters  listed  above 
can  be  computed  from  the  near-field  statistical  parameters. 
Thus,  pulsed  antennas  characterization  does  not  require 
knowledge  of  any  additional  statistical  near-field  data. 

(4)  The  probability  density  function  (p.d.f.)  for  the  far-field 
electric  fields  of  correlated  random  sources  or  scatterers 
can  be  accurately  approximated  for  most  antenna  or  scattering 
problems  via  the  use  of  Nakagami ' s  p.d.f.  and  linear  operator 
theory.  A  less  accurate  but  relatively  simple  expression 
for  the  p.d.f.  was  also  derived  that  is  suitable  for  certain 
engineering  applications. 

These  results  imply  that  the  fundamental  technical  requirement  for 
employing  near-field  techniques  to  describe  wideband  CW  or  pulsed 
antenna  radiation  is  that  one  be  able  to  compute  the  statistical 
average  far-field  power  pattern  versus  frequency  from  the  measured 
near-field  data.  As  stated  in  Item  2,  this  requires  a  knowledge  of  the 
average  electric  field  and  the  standard  deviation  at  each  sample  point 
and  the  covariance  function  at  all  different  measurement  points  (cross¬ 
covariance).  Only  the  cross-covariance  function  presents  a  significant 
measurement  problem.  The  accuracy  and  feasibility  of  the  near— field 
measurement  technique  for  wideband  out-of-band  antennas  depends  on 
whether  the  covariance  functions  can  be  suitably  determined. 
Accordingly,  considerable  effort  was  devoted  to  studying  the  covariance 
functions  and  their  effect  on  accuracy,  time  and  cost.  The  results  of 
this  effort  are  presented  and  discussed  in  subsection  C.  In 
particular,  the  theoretical  and  numerical  analyses  are  presented  for  a 
linear  array  of  wire  elements.  The  theory  and  equations  for  a  two 
dimensional  array  of  wire  elements  involve  no  new  concepts. 
Corresponding  theoretical  and  numerical  analyses  of  an  array  of 
multimoding  out-of-band  waveguide  elements  can  also  be  performed  by 
including  intermodal  correlations  in  the  analysis. 
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B.  Frequency  Domain  Analysis 


The  analysis  was  conducted  initially  for  the  linear  array  of 
waveguide  elements  depicted  in  Figure  2-1.  The  array  is  assumed  to  be 
either  pulsed  or  operated  oyer  a  finite-'  frequency  band.  The  frequency 

i 

spectrum  is  assumed  to  contain  some  components  which  can  be  propagated 
in  higher-order  mode  field  configurations  in  the  waveguide  elements.  A 
convenient  starting  point  in  the  analysis  is  to  first  write  down  the 
appropriate  equations  for  a  non-random  antenna  measurement  situation. 
The  analysis  of  randomly-excited  antennas  can  then  subsequently  be 
conducted  based  on  the  initially  deterministic  equations  and 
statistical  analysis  techniques. 

1.  Basic  Equations 

The  electric  field  produced  on  the  near-field  measurement 
plane  is  written  as  the  superposition  of  the  radiation  fields  of  the 
individual  element  radiators.  It  is  assumed  that  the  measurement  plane 
is  located  at  a  distance  x  that  is  ^  the  far-field  distance  of  each 
element  radiator.  A  transverse  field  component,  say  the  vertically 
polarized  component  Eq(<s),  produced  at  the  measurement  point  q  is  then 

K  K  ^  I 

E,«o)  .  ) - .  (2-1) 

£q 

where  c  is  the  speed  of  light  in  vacuum,  u  =  2TTf,  and  where 


=  complex  mode  coefficient  for  the  mode  in  the 
element, 

K 

h^  =  vertically  polarized  far-field  electric 
field  pattern  of  the  mode  of  the 
element, 


(J) 


i,q 


=  angular  location  of  the  q*-h  measurement  point 
with  respect  to  the  center  of  the  element. 


r.  ”  magnitude  of  the  radius  vector  from  the  center 
of  the  element  to  the  measurement 
point . 


The  plots  of  amplitude  and  phase  shown  near  the  bottom  of  Figure  1 
depict  either  (1)  the  electric  field  versus  frequency  as  the  array 
input  signal  is  swept  over  a  specified  frequency  band  or  (2)  the 
complex  frequency  spectrum  of  a  I'adiatcd  pulse.  In  either  case, 
different  results  would  be  obtained  at  different  measurement  points. 
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Figure  2-1.  Sketch  depicting  a  linear  array  of  waveguide  elraents  and 
hypothetical  out-of-band  aaiplitude  and  phase  responses 
at  a  near-field  aeasureaent  point. 


The  frequency  spectrum  of  the  time  pulse  is  obtained  via  the  Fourier 
transform  of  the  pulse,  which  is 


Eq  (o)) 


L  En(tn)  exp 
n  T 


(2-2) 


where  EqCt^)  is  the  time-domain  response. 

The  far-field  electric  field  is  obtained  as  the  discrete  Fourier 
transform  of  the  near-field  electric  field.  Thus,  the  far-field 
electric  field  is  given  as 


Z 

q 


Eq(w)  exp 


(2-3) 


where  Yq  is  the  y-coordinate  of  the  measurement  point,  (j)  is  the 

azimuth  angle  of  the  far-field  observation  point,  and  where  Eq  is  a 
previously  defined  by  Equation  (2-1).  A  factor  (l/r),  where  r  is  the 
distance  from  the  center  of  the  measurement  plane  to  the  far-field 
observation  point,  has  been  suppressed  in  Equation  (2-3)  and  subsequent 
equations . 

The  power  density  in  the  far-field  of  the  antenna  is  obtained  by 
multiplying  Equation  (2-3)  by  its  complex  conjugate.  The  resulting 
equation  for  the  power  density  P(w,(J>)  is  thence 


p((i),<t>)  =  E*(w,((>)E(a),<ti)  = 

^q'^q  E*q*  (w)Eq((D)  exp 


[j  ^  8in(')))  (Yq'-Yq)] 


(2- A) 


Equations  (2-1)  through  (2-4)  are  the  well-known  basic  equations 
for  analyzing  deterministic  (non-random)  antenna  patterns  utilizing  the 
frequency  domain  approach.  The  temporal  (time)  behavior  of  the 
electric  field  is  obtained  via  the  Fourier  transform  with  respect  to 
frequency.  The  corresponding  analysis  of  a  randomly-excited  array 
antenna  primarily  involves  the  application  of  certain  mathematical 
operations  to  these  same  equations,  as  described  in  the  following 
paragraphs. 

The  plots  shown  in  Figure  2-1  can  be  interpreted  as  representing 
the  measured  response  obtained  from  one  experiment  involving  a 


11 


r«ndo«il ^-excited  ar.tenna.  Suceeasive  experiaanta  would  yield  different 
responses.  Consequently,  a  randoaly  excited  antenna  has  aany  possible 
near-field  distributions,  spectral  responses,  and  far-field  patterns. 
Thus  the  anount  of  data  required  to  characterise  a  randoaly  excited 
antenna  will  be  such  greater  than  the  aaount  of  data  required  for  a 
comparable  deterministic  antenna  unless  a  suitable  statistical  analysis 
can  be  devised  to  reduce  the  data  requirements. 

The  statistical  average  value  of  the  far-field  power  density  is 
written  as 

<P(u,(f)>  - 


5:^,  <  E*q  <  (  w,  (>)  Eq(a),(>)>  exp  |^j  ^  sinO  )  (Yq<-Yq>j, 

(2-5) 

where  the  angular  brackets  denote  the  statistical  average  value  (  8l. 
The  angular  brackets  are  shorthand  notation  for  integrals  of  the  type 


<W(x)>  * 


-C  )d  j;  -- dc 
n  1  n 


(2-6) 


^1-^0 


where  x  ”  a  non-random  variable, 

Q  "  random  variables,  and 
f(4)  “  probability  density  for  the  random  variables. 

Similarly,  the  statistical  average  value  of  the  (complex)  electric 
field  is 


E(u,$)  ■  ^q<Eq(u>)>  exp  ^-j  ^  sin  (f)  Yqj  .  (2-7) 

The  electric  field  is,  of  course,  a  coi^lex  valued  function  and  is 
therefore  not  an  observable  quantity.  Bowever,  the  statistical  average 
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value  and  higher-order  statistical  moments  are  well  defined  quantities 
[  5  ].  In  particular,  the  statistical  average  value  of  the  complex 
electric  field  Eq  is  defined  as 

<Eq(u>)>  -  <Uq(u»)>  -  j  <Vq(u))>  .  (2-8) 

where  Uq  and  Vq  denote  the  real  and  imaginary  parts,  respectively,  of 
Eq.  Uq  and  Vq  are  defined  in  the  conventional  manner  as 

Uq(<^)  “  Aq(w)  cos  [Oq(w)]  ,  and  (2-9) 

Vq(‘*>)  -  Aq(a))  sin  (Oq(w)]  ,  (2-10) 

where  Aq(u))  •  relative  amplitude  (real  number)  of  the  electric 

field  at  q,  and 

o(q(‘^)  **  relative  phase  of  the  electric  field  at  q. 

Thus,  a  knowledge  of  the  average  values  of  the  real  and  imaginary  parts 
of  the  near-field  electric  field  over  the  measurement  plane  permits  the 
computation  of  the  average  far-field  electric  field  as  the  Fourier 
transform  of  the  average  near-field  electric  field. 

The  average  power  density  is  related  to  the  product  <  E*  >  <E>  , 

where  the  symbol  *  denotes  complex  conjugation,  as 

<P((D,<|>)>  -  <E*(w,$)><E(w,<{»)>  +  Ce*,e(**‘*'I*)  •  (2-11) 

The  symbol  Ce*,e  denotes  the  covariance  function  and  is  defined  as 

Ce*,e(“»'^>  "  <E*(<*',<fr)E('0,«fr)>  -  <E*(u,*)><E(w,'fr)>(2-12) 

The  equation  for  the  far-field  covariance  function  can  be  derived  for 
the  linear  array  via  straightforward  algebraic  manipulations.  The 
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resultlfig  equacioa  is 


(w.^>  + 

\  fj  7  sln{«)  <V^,-Y^)|  . 

Cq’Vq) 

(2-13) 


where  denotes  the  atsndard  deviation  of  the  near-field  electric 
field  at  q  and  Rq'q  denotes  the  cross-correlation  coefficient  for  the 
electric  field  at  q  and  the  conjugate  of  the  electric  field  at  q' . 

The  standard  deviation  Yq  is  a  real  number  and  is  equal  to  the 
square  root  of  the  sue  of  the  variances  of  the  real  and  imaginary  parts 
of  Eq.  Accordingly,  Yq  is  written  as 

Yq(w)  =  ^  ■*'  ^^q  ’  (2-lA) 


where  ^  standard  deviation  of  the  real  part  of  £q,  and 

}q  ~  standard  deviation  of  the  imaginary  part  of  Eq. 

The  cross-correlation  coefficients  Eq'q  are  defined  as  the  complex 
numbers  obtained  via  the  equation 


R  t  (‘^)  “  — 7— r 

q  q  Y^i  (<s)  r^(<»)) 


(2-15) 


where  the  numerator  of  Equation  (2-15)  is  referred  to  as  the  cross- 
covariance  function.  The  croas-covariaoce  of  the  electric  fields  at  q 

and  q'  is 
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(2-16) 


C  ,  (a.)  =  j 

q  q  1 

J<U^,(ai)U^(u))>  -  <U^,(u>Xu^(w)>j 

■h 

f<Vq,(w)V^(m)>  -  <V^.(w)Xfq(w)>j| 

-1 

j<V^,(a))U^(m)>  -  <V^,(uiXJ^(a))>j 

- 

|<U^,(a))V^(m)>  -  <tJq,('*)><V^(a)t>J  j 

The  rigorous  analysis  of  the  far-field  statistical  average  power 
pattern  for  a  given  frequency  is  seen  from  Equation  (2-7)  through  (2- 
16)  to  require  computation  of  the  following  near-field  statistical 
quantities : 

(1)  statistical  average  value  of  the  real  and  imaginary  parts 
of  the  near-field  electric  field  at  all  measurement  points, 

(2)  the  standard  deviation  of  the  real  and  imaginary  parts  of 
the  near-field  electric  field  at  all  measurement  points,  and 

(3)  the  cross-polarization  coefficients  of  the  near-field 
electric  field  at  all  different  measurement  points. 

In  a  near-field  measurement  situation,  the  first  two  quantities  can  be 
determined  by  computing  the  "sample"  average  values  and  standard 
deviations  obtained  from  repeated  trials.  The  sample  average  value  [8  ] 
of  a  random  variable  W  is  defined  as 

<">■  I  '^n"n  • 

where  is  the  value  of  W  obtained  in  the  n^^  trial,  and  N  is  the 
number  of  trials.  Similarly,  the  sample  standard  deviation  ow 
defined  as 
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(2-18) 


3  « 

w  “ 


[H^  -  <W>r 


The  extraction  of  the  cross  correlation  data  is  not  straightforward  and 
this  probl.ea  area  is  discussed  in  Part  C. 

2.  Frequency-Averaged  Statistical  Average  Antenna  Patterns 

The  preceding  discussions  have  addressed  statistical 
averaging  over  randomly -varying  variables.  It  may  be  meaningful  to 
also  average  over  frequency  in  some  applications  involving  a  CW 
radiating  system  of  aioderate  bandwidth.  The  purpose  of  averaging  over 
frequency  is  to  obtain  a  single  average  pattern  plus  standard  deviation 
that  adequately  describes  the  general  radiation  characteristics  of  Che 
antenna.  The  single  average  pattern  plus  standard  deviation  replaces 
the  large  collection  of  patterns  versus  frequency  that  would  otherwise 
be  needed  to  characterize  Che  antenna. 

Two  different  methods  for  obtaining  the  frequency-averaged 
statistical  pattern  have  been  formulated.  The  most  direct  method  for 
obtaining  the  frequency-averaged  statistical  average  pattern  is  to 
first  compute  the  statistical  average  pattern  at  selected  frequencies 
and  then  arithmetically  average  the  statistical  average  patterns.  This 
process  is  described  mathematically  as 

=  I  .  (2-19) 

where  N  is  the  total  number  of  selected  frequencies  and  where  frequency 
averaging  is  denoted  by  subscript  u>  on  the  outermost  right-hand 
angular  bracket.  It  is  also  possible  to  obtain  the  frequency-averaged 
statistical  average  pattern  by  first  computing  the  deterministic 
frequency  average  and  then  computing  the  statistical  average.  This 
process  is  described  mathematically  as 

<P(u). <!>)>>  -  r.  ,  L  <D  ,  >  exp[-jd  0)].  C2-20'J 


If. 


where 


d  .  -  (Y  -  Y  ).  and 

q  q  c  q  q 


(2-21) 


(2-22) 


The  aaount  of  conputational  labor  required  to  coapute  frequency- 
avorage  statistical  average  patterns  is  roughly  equivalent  for  the  two 
■ethods. 

The  out-of-band  radiation  patterns  of  a  20-eleaent  BwltiBoding 
array  were  studied  numerically  via  both  direct  Monte  Carlo  simulations 
and  Equations  (2-1)  through  (2-19).  For  this  array,  whose  in-band 
design  frequency  is  9.0  GHz,  Monte  Carlo  patterns  as  well  as  the 
statistical  average  patterns  and  standard  deviations  were  computed  for 
selected  frequencies  up  to  19  GHz  for  various  in-band  scan  conditions. 
Ezperistentally-derived  out-of-band  phase  shift  and  attenuation 
statistical  parameters  for  a  waveguide  element  containing  a  ferrite 
phase  shifter  were  used  as  inputs.  Accordingly,  the  modal  phase 
variations  follow  a  Gsussian  distribution  and  the  modal  power 
variations  follow  a  unifcrm  distribution  to  good  approximation. 

It  was  found  that  for  small  in-band  scan  angles  the  statistical 
average  pattern  and  standard  deviation  at  frequencies  of  14  GHz,  15.7 
GHz,  and  18  GHz  could  provide  a  good  EMC  description  of  the  radiation 
characteristics  over  portions  of  the  frequency  band  for  which  energy 
can  propagate  in  the  first  two  modes,  the  first  three  modes,  and  the 
first  five  modes,  respectively.  Figure  2-2  shows  the  calculated 
average  pattern  and  standard  deviation  at  18  GHz  for  broadside  scan 
superimposed  on  calculated  Monte  Carlo  patterns  for  16-19  GBz  region 
for  the  cross  polarized  component  of  the  radiated  power  for  propagation 
in  the  first  5  waveguide  sndes.  For  in-band  scan  angles  of  about  30 
degrees  or  more,  a  better  EMC  description  is  obtained  by  including  at 
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least  the  statistical  average  patterns  for  frequencies  near  the  edges 
of  the  frequency  sub-interval  because  the  patterns  at  different  out-of- 
band  frequencies  scan  through  unequal  amounts  for  a  given  in-band  scan 
angle  [  9  ].  The  spatial  regions  corresponding  to  the  superimposed  main 
beams  and  grating  lobes  clearly  delineate  areas  of  particular  concern 
in  EMC  applications. 

The  statistical  average  patterns  and  standard  deviations  at 
selected  frequencies  can  provide  a  very  succinct  engineering 
description  of  the  important  EMC  characteristics  of  wideband  CW 
multimoding  antennas.  The  statistical  average  patterns  and  standard 
deviations  are  a  practial  alternative  to  the  comparatively  more 
expensive  and  cumbersome  Monte  Carlo  simulations. 

3.  Reduction  of  Data  Measurement  Requirements 

A  reduction  in  the  near-field  data  measurement  requirements 
can  be  achieved  at  the  expense  of  a  nominal  reduction  in  the  accuracy 
of  the  far-field  pattern  details.  Specifically,  reducing  the  amount  of 
measured  data  either  by  increasing  sample  spacing  or  by  truncating  the 
size  of  the  transverse  measurement  plane  will  generally  cause  errors  in 
the  calculated  far-field  patterns.  The  magnitude  of  the  far-field 
pattern  errors  is  a  function  of  the  sample  spacing,  the  relative  power 
level  at  the  edges  of  the  truncated  near-field  measurement  plane,  and 
the  accuracy  of  the  measured  amplitude  and  phase  recorded  at  each 
sample  point  [lOl.  The  effects  of  number  of  sample  points  and  phase 
measurement  accuracy  on  the  calculated  far-field  power  pattern  levels 
is  displayed  in  Figure  2-3.  The  effect  of  the  number  of  sample  points 
on  the  beam  pointing  error  are  displayed  in  Figure  2-4  for  a  specified 
phase  measurement  error.  It  was  assumed  that  the  relative  power  at  the 
edges  of  the  measurement  plane  is  -25  dB  with  respect  to  the  highest 
recorded  near-field  amplitude.  The  analysis  and  equations  employed  to 
obtain  these  results  are  summarized  in  Reference  [lOl. 

Figure  2-3  is  a  plot  of  the  error  in  the  calculated  far-field 
power  pattern  versus  the  "no-error"  power  pattern  level  for  an  R.M.S. 
have  measurement  error  of  5.0  degree  for  sample  spacings  of  IX  and  2^, 
where  ^  is  the  operating  wavelength.  The  "no-error"  power  pattern 
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RMS  BEAM  POINTING  ERROR  (DEGREES 


RMS  PHASE  ERROR  (DEGREES) 


Figure  2-4.  RMS  beaa  pointing  error  as  a  function  of  RMS 

phase  error  for  an  antenna  whose  3-dB  beaawldth 
is  noainally  0.5  degrees  for  the  indicated  sample 
spacings. 


level  along  the  abscissa  is  the  power  pattern  level  calculated  from  X/2 
spacing  for  a  near-field  measurement  plane  which  extends  to  the  -40  dB 
level.  Figure  2-4  is  a  plot  of  pointing  error  versus  R.M.S.  phase 
measurement  error  for  sample  spacings  of  X/2',  IX  ,  and  2X 

Clearly,  the  accuracy  of  the  calculated  power  pattern  is  degraded 
by  large  sample  spacings,  truncated  measurement  planes,  and  near-field 
measurement  errors.  It  is  also  true  that  smaller  near-field 
measurement  errors  permit  the  use  of  larger  sample  spacings  and/or  more 
severely  truncated  measurement  planes.  Of  course,  the  advantages  of 
reduced  data  recording  requirements  must  be  weighed  against  the 
accuracy  requirements  for  the  antenna  under  test.  However,  it  appears 
that  substantial  reductions  can  be  achieved  for  many  engineering 
applications  where  the.  main  interest  is  in  obtaining  valid  estimates  of 
mainbeam  and  grating  lobe  power  levels  and  pointing  directions  and 
where  moderate  over-estimates  of  the  average  sidelobe  level  are 
acceptable . 

C.  Near-Field  Covariance  Study 

The  fact  that  the  electric  fields  at  different  near-field  sample 
points  are  correlated  has  an  effect  on  the  measurement  time,  accuracy 
and  cost  which  can  be  appreciated  by  considering  two  different 
measurements  using  a  single  probe  to  measure  a  wideband  phased  array 
antenna.  The  two  measurement'  schemes  will  be  denoted  for  convenience 
as  Method  I  and  Method  II. 

In  Method  I,  data  are  recorded  as  the  probe  is  moved  over  the 
measurement  plane  in  a  raster  scan.  The  operating  frequency  and  array 
scan  condition  are  assumed  to  be  constant  while  the  probe  r^cans  the 
entire  measurement  plane.  This  takes  about  two  hours  of  elapsed  time. 
The  probe  is  then  returned  to  its  starting  point,  and  the  array  phase 
shifters  are  "cycled"  and  returned  to  the  initial  scan  condition.  Data 
are  recorded  at  the  same  frequency  as  the  probe  again  traverses  the 
measurement  plane.  This  sequence  of  events  is  repeated,  say,  50  times 
for  a  given  frequency  and  array  scan  condition.  The  total  measurement 
time  required  to  record  data  for  100  frequencies  and  30  array  scan 
conditions  is  about  300,000  hours.  The  measurement  time  is  clearly 
excessive.  Additionally,  it  is  very  unlikely  that  an  out-of-band  array 


22 


antenna  under  test  can  be  kept  stable,  i.e.,  no  variation  in  modal 
content  in  any  element,  for  the  two-hours  needed  to  record  data  for 
each  selected  out-of-band  frequency  and  array  scan  condition.  However, 
if  the  array  element  signals  could  be  kept  stable  the  near-field 
covariance  functions  could  be  computed  directly  from  the  measured  data. 

In  Method  II,  the  probe  is  held  stationary  at  a  preselected  sample 
point  while  the  signal  source  is  scanned  through  the  entire  frequency 
spectrum  of  interest  for  a  given  array  scan  condition.  The  phase 
shifters  are  then  cycled  and  returned  to  the  selected  scan  condition. 
This  procedure  is  repeated  50  times.  A  pew  scan  condition  is  selected 
and  the  foregoing  sequence  of  events  is  repeat-id.  Finally,  the  entire 
sequence  of  events  is  repeated  for  all  of  the  preselected  sample 
measurement  points  on  the  near-field  plane.  Total  measurement  time  for 
100  frequencies  and  30  array  scan  conditions  is  estimated  to  be  about 
120  hours.  This  measurement  time  is  reasonable  for  a  thorough 
characterization  of  a  randomly-excited,  out-of-band,  wide-band  phased 
array  antenna.  The  average  value  and  standard  deviation  of  the 
electric  field  at  each  near-field  sample  point  can  be  computed  directly 
from  the  recorded  data  for  each  frequency  and  scan  condition.  However, 
Che  covariance  function  for  the  electric  fields  at  different  near-field 
sample  points  is  more  difficult  to  extract  from  the  recorded  data. 

The  covariance  functions  for  the  near-field  electric  field  of  a 
given  array  cap  be  computed  in  a  straightforward  manner  in  terms  of  the 
covariance  function  for  the  electric  fields  on  the  array  aperture.  Of 
course,  the  aperture  covariance  function  is  not  known  in  practice.  If 
it  were  known,  there  would  be  no  need  for  near-field  measurements. 
However,  an  analysis  of  a  specified  array  antenna  can  be  used  to  study 
the  general  nature  of  the  near-field  covariance  functions  and  their 
impact  on  the  accuracy  of  the  far-field  average  power  pattern 
calculatiop.Hi .  Furthermore,  near-field  covariance  functions  computed  in 
this  manner  provide  baseline  data  that  can  be  used  to  check  the 
validity  and  accuracy  of  near-field  covariance  functions  obtained  from 
analysis  of  simulated  measured  data  obtained  from  a  simulated  Method  II 
measurement.  Accordingly,  the  near-fieid  covariance  functions  for  a 
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1  inear  array  of  wire  dipoles  have  been  studied  analytically  and 
numerically.  A  corresponding  study  for  an  array  of  multimoding 
waveguide  elements  would  be  beneficial  but  is  beyond  the  scope  of  the 
current  efforts. 

1.  Analysis  of  Wire  Array 

The  linear  array  of  vertically-oriented,  center-fed  wire 
elements  shown  in  Figure  2-5  is  analyzed  in  the  following  paragraphs. 
It  is  assumed  in  the  analysis  that  each  element  is  fed  by  a  constant 
voltage  source. 

The  electric  field  Eq  at  measurement  point  q  on  the  line  defined 
by  Z=0,  X=Xo  is  Z-directed  and  may  be  written  as 


E 

q 


c  X 

o 


V^o  * 


(2-23) 


where 

=  complex  current  at  the  current  amplitude  maximum 
on  the  i  element, 

Yq  =>  the  Y  coordinate  of  the  measurement  sample  point, 
Yjj^  =  the  Y  coordinate  of  the  H  element, 
k  =  Itt/ X  where  X  is  the  wavelength,  and 


The  integral  appearing  in  the  definition  of  Cq  is  the  same  for  each 
element  at  a  given  frequency  and  is  just  the  Fourier  Transform  of  the 
current  distribution  (f(z')  along  the  length  of  the  element,  where 
^(z')  is  normalized  to  have  a  maximum  amplitude  of  1.0.  Cq  will  be 
suppressed  in  most  of  the  remaining  equations. 

The  statistical  average  far-fiold  power  density  ?(<{>)  at 
azimuth  angle  <{>  can  be  computed  as  previously  stated,  to  wit: 
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NEAR-FIELD 
SAMPLE  POWTS 


Figure  2~5.  Sketch  of  an  array  of  center-fed  vire  elenents  and 
near-field  aaa^le  points. 
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<P($)>  -  <£*(♦)>  <E(#)>  ♦  exp  [(jk  »in  #)  (Yq*  -Yq)l 

(2-24) 


where  <£(#)>  is  the  Fourier  Trsnsfons  of  the  sverage  near-field 
electric  field,  <^E*(^)I>  is  its  coaplex  conjugate,  and  Cq>q  is  the 
near-field  covariance  function  defined  as 

Cq'q  *  <Pq*  ®q>  '  <E^’>  <Eq>  (2-25) 

%>here  the  sysibol  *  denotes  coaplex  conjugation. 

The  near-field  covariance  function  Cq'q  obtained  froa  Equations 
(2-23)  and  (2-25)  is  then  explicitly  coiiputed  as 


t'  t 


1*1  I'q*  Iq 


(2-26) 


where 

C  ■  the  covarance  function  for  the  current  on 
eleaent  1  and  the  conjugate  of  the  current 
on  eleaeot  1 ' ,  and 


The  current  covariance  fu.tction  G  is  by  definition 

‘^I’l  “  ^1^  ■  •  (2-28) 

Equations  (2-26)  and  (2-28)  express  Che  fact  Chat  Che  near-field 
covariance  function  can  be  coBq>uted  froa  Che  current  covariance 
function. 
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The  current  covariance  function  is  a  function  of  the  oiutual 

ii'hnUtatite  Matt  L>  for  the  array  anti  the  atatlstiral  paraoiateta  of  the 
input  voltage  sources.  In  particular,  Che  current  aC  the  element 

is  computed  as 


^  r  V 

£  K  £k  K 


(2-29) 


where 


elements  of  Che  complex  admittance  matrix  [11],  and 
V  “  complex  voltage  of  the  ic^h  voltage  source. 


It  follows  Chen  Chat  C ^  is  computed  as 


c „ , „  E  , E  r„ ,  ,r,  c  , 

£*£  ic'k£'k*£k:  tc'i 


(2-10) 


where  is  Che  voltage  covariance  function. 


The  current  covariance  function  has  zero  magnitude  in  two 
different  special  cases.  First,  the  current  covariance  function  is 
zero  when  all  of  the  voltage  covariance  functions  are  zero.  This 
occurs  when  the  voltages  are  deterministic  (non-random).  The  currents 
are  then  perfectly  correlated  with  correlation  coefficient  R£,£  *  1*0 
as  can  be  discerned  from  the  equation 


£'£  C*C£ 


(2-31) 


where  *  standard  deviation  of  the  complex  current,  by  taking  the 

limit  as  the  numerator  and  denominator  approach  zero.  Second,  the 
current  covarionce  function  for  Z  ^  Z'  is  zero  when,  simultaneously, 
the  voltages  ate  uncorrelated  and  the  mutual  coupling  between  current 
elements  is  zero.  In  this  case  the  voltage  covariance  C  ,  **  6  , 

X  K  K  X 

(Yy)|^  and  ^£x  -  the  voltage  covariance  matrix  and  the 
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mutual  admittance  matrix  both  contain  non-z^ero  matrix  elements  only 
along  the  diagonal.  This  second  case,  where  C  *  0  for  t  n' , 

Jv  X 

corresponds  to  the  situation  where  the  currents  on  different  elements 

are  uncorrelated  and  hence  C  .  “  6  .  . 

£*i  1 

It  will  be  assumed  in  the  rest  of  the  analysis  that  the  voltage 

sources  are  statistically  independent  and  are  therefore  uncorrelated. 

The  behavior  of  the  near-field  covariance  function  will  be  examined  for 

the  case  involving  isolated  current  elements  (r.  *=  6  T  )  and  the  case 

ttc  lie  1 

involving  mutual  coupling  among  the  current  elements. 

The  near-field  covariance  function  may  be  written  as 


Cq'q 


1  11' 


^^v^l  ®lq’  ^Iq 


^‘^1  '/i-/ik<Vk 


“^I'lC^lK^  ^I'q'^lq  (2-32) 


for  the  realistic  situation  when  the  currentc  are  coupled  and  as 


-q-q 


^■'"v^l  *'lq’®lq 


(2-33  ) 


when  the  element  currents  are  not  coupled.  Equation  (2-33)  is  obtained 

from  Equation  (2-32)  by  setting  r**  T.  "6.6  so  that  the  o££- 

l  K  IK  I’k  Ik 

diagonal  elements  of  the  mutua'  -*dmittance  matrix  are  zero.  Equation 
(2-33)  shows  that  the  near-fie  covariance  function  is  non-zero  even 
if  the  near-field  electric  field  is  produced  by  isolated,  statistically 
independent  currents.  This  covariance  will  be  denoted  for  convenience 
as  intrinsic  covariance.  It  is  always  present  in  stochastic  antenna 
problems.  Equation  (2-32)  shows  that  the  total  near-field  covariance 
function  in  the  realistic  situation  involving  coupled  radiating 
elements  consists  of  two  components,  namely  the  intrinsic  covarisnee 
and  a  second  component  denoted  herein  as  the  interactive  covariance. 
The  interactive  covariance  is  present  only  sdieo  the  radiating  elements 
"interact",  i.e.,  are  electromagnetically  coupled.  A  third  component 
of  the  nesr-field  covariance  function,  denoted  as  the  intra-element 
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covariance,  occurs  for  aulciaoding  waveguide  antennas  due  to 
correlation  anong  the  nodes  propagating  within  each  separate  waveguide 
eleiaeat  proper.  This  covariance  contribution  is  absent  for  the  wire 
array  under  conaideratioc. 

It  should  be  noted  that  the  near-field  covariance  function  for  the 
wire  array  given  by  Equation  (2-32)  is  scan- invar iant .  This  result 
follows  froB  the  fact  Chat  the  nutual  adnittances  are  scan-invariaut . 
Of  course,  the  statistical  average  asq^litudes  and  relative  phases  of 
the  electric  field  at  each  sanple  point  do  change  with  scan  angle.  The 
near-field  covariance  functions  for  a  nultinoding  out-of-band  y>avegulde 
array  nay  change  with  scan  angle  due  to  changes  in  average  node 
excitations  with  scan  angle. 

2.  Nunerical  Study  of  Wire  Array 

Nuaierical  ainulations  were  conducted  to  study  both  the  near- 
field  and  far-field  radiation  charactei istics  of  a  linear  array  of  nine 
center-fed  wire  eleaents.  The  array  elenents  have  length  Aq/2,  where 
is  the  free-apace  wavelength  at  the  design  in-band  frequency  of  3.0 
GBz,  and  are  spaced  1^/2  apart  along  Che  y  axis.  Each  wire  elenenC  is 
assumed  to  be  fed  by  a  constant  anplitude  voltage  source,  and  each 
source  is  assusMid  to  be  sutched  to  a  50-(Hib  at  both  in-band  and  out-of- 
band  frequencies. 

The  following  data  were  cosputed  for  an  array  of  coupled  elenents 
and  an  array  of  "isolated"  eleaents  for  both  the  in-band  frequency  of 
3.0  GHz  and  and  the  out-of-band  frequency  of  9.0  GHz: 

(1)  Non-randon  (deteminisCic)  near-field  power  diatributionc 
and  far-field  power  paccems, 

(2)  Sandon  (Monte  Carlo)  near-field  power  distributions  and 
far-field  power  paCCerua, 

(3)  Analytical  staCiaCical  average  near-field  power 
distribution  and  analytical  staCiatical  average  far- 
field  power  pattern, 

(4)  NusMrical  averse  far-field  pattern  and  its  associated 
standard  deviation  based  on  50  McnCe  Carlo  randos 
far-field  patCeros,  and 

(5)  Hear-field  covarianca  functions. 


The  data  for  Itese  (1)  through  (A)  were  coaputed  for  in-band  acan- 
angles  of  zero  degrees  and  30  degrees.  The  near-field  covariance 
functions  for  the  wire  array  are  scan-invariant  and,  hence,  they  were 
explicitly  computed  only  for  the  broadside  (zero-degree)  scan 
condition. 

The  near-field  covariance  functions  and  their  effects  on  the 
statistical  average  power  patterns  are  the  topics  of  paramount  interest 
and  only  there  data  will  be  presented  and  discussed  herein.  A  sketch 
of  the  simulated  near-field  sieasurement  situation  is  shown  in  figure 
2-6.  The  near-field  data  were  computed  for  65  saBq>le  points  centered 
on  the  near-field  sample  line.  The  near-field  sample  line  was  chosen 
to  be  8 A  wide  and  is  located  2  A  from  the  array. 

The  key  results  of  the  numerical  studies  can  be  auamarized  with 
the  aid  of  the  data  plotted  in  the  group  of  Figures  2-7  through  2-10 
for  the  in-band  frequency  of  3.0  GHz  and  the  group  of  Figures  2-11 
through  2-lA  for  the  out-of-band  frequency  of  9.0  GHz.  Thi  data  are 
sequenced  in  the  same  order  within  each  group.  The  first  figure  in 
each  group  is  a  plot  of  the  near-field  power  distribution  for 
deterministic  (non-random)  excitation  of  an  array  of  coupled  current 
eletsents  for  broadside  scan.  This  plot  is  followed  by  plots  of  the 
scan-invariant,  near-field  statistical  covariance  functions  and  plots 
of  the  corresponding  statistical  average  far-field  power  patterns  for 
broadside  scan.  The  statistical  data  In  each  figure  are  displayed  for 
interacting  (coupled)  array  elements  and  non-anteracting  (isolated) 
array  elements  by  the  plots  labeled  B  and  A,  respectively.  All  of  the 
nesr-field  plots  in  Figures  2-7  through  2-lA  are  normalized  with 
respe-t  to  the  peak  magnitude  of  the  non-random,  in-band,  near-field 
power  distribution  shorn  in  Figure  2-7. 

The  general  behavior  of  the  near-field  covariance  functions  may  be 
discerned  from  Figures  2-8  and  2-9  for  the  in-band  frequency  and 
Figures  2-12  and  2-13  for  the  ouC-of-band  frequency.  The  covariance 
functions  for  the  center  sample  point  with  all  other  points  are 
symmetric  about  the  center  point,  while  the  covariance  funccicns  for  a 
sample  point  opposite  the  edge  of  the  array  with  all  other  saq>le 
points  are  asymmetrical.  This  behavior  is  expected  because  the  arrey 
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Figure  2-7.  Non-random  near-fleld  power  di8tribut<on  for  the  In-baod  frequency 
of  3.0  GHz  for  the  In-band  scan  angle  of  0  degrees. 


Figure  2-8.  Intrinsic  (A)  snd  total  (fl)  near-£ield  covariance  function  of  the 
eldccrlc  field  at  y*0  and  the  conjugate  of  the  electric  field 
at  all  other  points  on  the  near-fleld  plane  for  the  in-band 
frequency  of  3  GHs  for  the  In-band  scan  angle  of  0  degrees. 


«r  all  ocher  points  on  the  neor-fleld  plane  for  the  In-band 
frequency  of  3  CK*  for  the  In-band  scan  aaifle  cf  0  degrees. 


Figure  2-10.  Analytical  atatlstical  average  pattern  for  a  nlne-elenent  dipole 
array  of  interacting  <B)  and  non-interacting  (A)  eleaents  for 
the  iit'-baad  frequency  of  3.0  GHz  for  the  in-hand  scan  angle  of 
0  degreea. 


rigure  2-11.  Kon-r*ndo«  near-fleld  power  distribution  for  the  out-of-band  frequency 
of  9.0  GHz  for  the  in-band  scan  angle  of  0  degrees. 
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Figure  2-13.  Intrinsic  (A)  and  total  (B)  near-fleld  covariance  function  of  the 
electric  field  at  y--20  and  the  conjugate  of  the  electric  field 
at  all  other  points  on  the  near-fleld  plane  for  the  out-of-band 
frequency  of  9  GHz  for  the  in-band  scan  angle  of  0  degrees. 
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is  finite.  An  infinite  array  would  everywhere  have  a  syanetriccal 
covariance  function.  The  total  covariance  function  has  greater  peak 
nuignitude  chan  Che  intrinsic  covariance,  and  also  has  broader  lobes. 
These  differences  are  more  pronounced  for  the  in'-band  frequency.  This 
is  reasonable  because  the  current  elesMnCs  are  only  weakly  coupled  at  9 
GHz.  Accordingly,  the  interactive  covariance  component  of  the  total 
covariance  is  much  smaller  Chan  Che  intrinsic  covariance  component  anu 
hence  Che  total  and  intrinsic  covariances  do  not  differ  markedly  at  9.0 
GHz. 

Ihe  behavior  of  the  statistical  average  far-field  power  patterns 
may  be  examined  vii  Che  plots  shown  in  Figure  2-10  for  3.0  GHz  and 
Figure  2-lA  for  9.0  GHz.  The  statistical  average  patterns  for  the 
array  of  interacting  current  elements  have  greater  peak  magnitu^'e  and 
decrease  more  rapidly  versus  azimuth  angle  than  the  average  patterns 
for  Che  array  of  non-interacting  current  elements.  The  difference  in 
Che  average  partcerns  for  interacting  and  non-interacting  elements  is 
more  pronounced  for  the  in-band  frequency  than  for  the  out-of-band 
frequency.  This  result  could  have  been  anticipated  from  the  near-field 
covariance  functions  because  Che  total  and  intrinsic  near-field 
covariances  are  significantly  different  for  3.0  GHz  but  are  only 
slightly  different  for  9.0  GHz.  The  results  for  other  scan  angles  not 
displayed  herein  are  entirely  consistent  with  the  results  for  the 
broadside  scan  angle. 

It  can  be  inferred  from  these  results  Chat  the  intrinsic 
covariance  function  is  a  "good"  engineering  approximation  to  the  total 
covariance  function  if  the  inter-element  coupling  is  sufficiently  weak. 
Further  numerical  studies  or,  preferably,  measured  data  are  needed  to 
determine  if  the  inter-element  coupling  of  multimoding  out-of-band 
waveguide  phased  arrays  is  sufficiently  weak  to  permit  valid  eaCimmtea 
of  the  average  far-field  power  patterns  to  be  obtained  in  this  manner. 
Of  course,  it  is  preferable  to  determine  the  total  covariance  functions 
from  Che  measured  data.  This  possibility  is  addressed  in  the  following 
subsection . 
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3.  Recovery  of  Covariance  Puoctioot  from  Measured  Data 

In  principle,  the  cover ience  functions  cen  be  recovered  froa 
the  Beeaured  steatieticel  data  via  Matrix  algebra,  aa  followa.  The 
near-field  covariance  function  ts  related  to  the  current 
covariance  function  via  Equation  (2-26),  to  wit: 


'q  q 


t’q'  °tq 


(2-34) 


where  G.  and  G  ,  ,  are  the  Green's  function  and  its  conjugate,  as 
previously  defined  by  Equation  (2-27).  This  equation  is  equivalent  to 
the  matrix  equation 


AD  -  B 


(2-35) 


where  D  and  B  are  column  vectors  corresponding  to  the  eleaent  current 
covariances  and  the  near-field  covariances,  respectively,  and  where  A 
is  the  matrix  formed  from  the  product  of  the  Green  function  and  the 
appropriate  conjugate  Green  function.  Equation  (2-35)  may  be  written 
in  terms  of  the  autrix  eleewnts  A^j  and  the  column  vector  components  Dj 
and  8^  as 


Ej  Aij  Dj  -  Bi  (2-36) 

where  Che  indices  i  and  j  are  related  to  the  (i',£)  indices  and  (q',q} 
indices  as 


i  ■  q  • 


(2-37) 


where  q'  is  fixed,  and 


j  •  t  ♦  (£•  -  1)  L. 


(2-38) 
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If  A  is  square  and  non-singular  and  the  right  hand  coluam  vector  B 
is  known  froa  aeasuresunt,  the  unique  solution  for  2>  can  be  obtained  by 
solving  the  set  of  linear  equations  obtained  frow  Equation  (2-36).  The 
index  j  aust  vary  from  1  to  where  L  is  the  nuaber  of  radiating 

elements,  and  ve  choose  i^l  to  so  that  A  will  be  a  square  aatrix. 
This  choice  cf  oaxiaum  i'^L^  aeans  that  aeasured  values  of  the  selected 
near-field  covariance  function  aust  be  supplied  for  saaple  points  on 
the  near-field  plane. 

The  matrix  A  will  be  non-singular  if  its  determinant  is  non-zero 
[12].  Computer  calculations  of  the  determinants  for  B-eleaent  array 
and  for  a  9-eleaenc  array  were  perfonaed  using  double  precision 

airthaetic.  These  coaputations  wore  conducted  for  aeveral  different 
ncar-field  distances  and  for  several  different  sets  of  saaple  point 
spacings  for  each  rtear-field  distance.  It  was  found  that  the  computed 
determinants  tor  the  S-eleaent  array  ranged  from  10~5  to  10~®.  The 
computed  determinants  for  the  9-eleaent  array  ranged  from  10~)2  to 
These  results  imply  that  the  matrix  A  is  "nearly"  singular. 
These  numerical  results  do  not  prove  chat  the  matrix  is  non-singular 
because  Che  compucacions  are  subject  to  non-negligible  errors  despite 

Che  use  of  double  precision  arithmetic.  However,  inspection  of  Che 

matrix  elements  further  supports  the  notion  that  the  matrix  is  non- 
singular  but  that  it  is  nearly  singular.  The  elements  of  Che  9x9 

matrix  for  the  3-elemenc  array  are  shown  in  Tsble  2-1.  It  can  be 

discerned  Chat  the  matrix  elements  do  not  dif'  markedly  from  each 
other.  However,  inspection  of  the  full  matriT  ows  that  no  two  rows 
are  Che  same  nor  are  any  two  rows  related  by  a  scant.  Similarly,  no 

two  columns  are  the  same  nor  are  any  two  columnc  ti  ’ted  by  a  constant. 
Thus,  one  may  conclude  that  the  matrix  is  non-singu.ar .  One  may  also 
conclude  from  the  determinant  computations  and  from  the  fact  Chat  the 
siatrix  elements  are  so  similar  that  the  matrix  is  "nearly"  singular. 

One  anticipates  trouble  in  numerical  computations  involving  a 
nearly  singular  matrix.  Most  computer  programa  designed  Co  numerically 
perform  utrix  inversions  or  to  solve  systems  of  equations  work  well 
for  a  diagonally  dominant  matrix,  but  they  can  fail  badly  for  a  matrix 


TKE  9X9  HATRIX  ARISING  IN  THE  SOLUTION  CP  THE 
NEAR-FIELO  COVARIANCE  FUNCTIONS  FOR  A  WIRE  ARRAY 
CONSISTING  OP  3  RADIATING  ELEMENTS* 
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that  departs  "too  much"  from  the  diagonally  dominant  situation.  The 
matrix  A  is  certainly  "ill  conditioned"  in  comparison  to  a  diagonally 
dominant  matrix.  Consequently,  the  more  common  methods  such  as 
Gaussian  elimination  with  iterative  improvement  or  the  iterative  Gauss- 
Seidel  and  Jacobi  techniques  [13]  failed  to  converge  to  accurate,  or 
even  sensible,  answers  for  the  9-element  array.  A  straightforward  row 
reduction  technique  written  for  this  problem  gave  good  results  for  3- 
element  and  5-element  arrays,  but  produced  grossly  erroneous  results 
for  7-element  and  9-element  arrays. 

Acceptable  numerical  results  were  finally  obtained  with  an 
algorithm  based  on  a  projection  method  for  solving  singular  or  nearly- 
singular  systems  of  equations  [  14  1.  This  same  method  has  been  employed 
in  the  analysis  of  probe  compensation  errors  conducted  under  this  same 
contract  and  is  well  documented  in  Reference  [15].  The  basic  method  of 
solution  can  be  summarized  by  considering  the  expression  for  the 
(m+l)th  iterate. 


D (m+ 1 ) 


D“  - 


m 


a_] 

m 


a 

m 


(2-39) 


where  d(°’)  is  the  solution  after  m  iterations  of  this  equation,  ag|  is 
the  m^^  column  vector  of  the  transpose  conjugate  matrix  of  A,  Bq  is  the 
mth  element  of  B,  and  where  the  square  brackets  denote  the  inner 
product  of  the  enclosed  quantities.  A  formal  proof  that  this  sequence 
converges  for  any  A,  B,  and  initial  guess  d(o)  is  presented  in 
Reference  [141.  If  the  system  of  equations  is  consistent,  it  then 
follows  that  the  limit  point  is  a  solution  of  the  system. 

The  real  parts  of  the  exact  values  of  the  normalized  current 
covariances  and  the  values  computed  from  the  projection  algorithm  are 
plotted  in  Figure  2-15  for  a  9-element  wire  array  operating  at  3.0  GHz. 
The  plotted  data  are  normalized  to  the  largest  value  of  the  exact 
current  covariance.  The  data  are  symmetrical  about  the  center  point 
(5,5)  and,  hence,  only  the  data  for  (1,1)  through  (5,5)  are  plotted. 
The  magnitudes  of  the  imaginary'  parts  of  the  exact  and  estimated 
covariances  are  0.012  and  0.04,  respectively,  and  are  therefore 
not  plotted.  The  computed  values  plotted  in  Figure  2-15  were  obtained 
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A  EXACT 


The  ex&ct  and  the  approximate  ulencnt  current  covariance  func 
for  the  9-eleQient  wire  array  lor  the  frequency  of  3-0  GHz  for 
uncorrelaced  Input  voltages. 


by  iterating  Equation  (2-39)  100  times,  which  took  57  seconds  of  cpu 
time.  Additional  iterations  produce  neglible  improvement  of  the  data. 
The  initial  starting  value  for  the  iterations  was 

=(L+/l  )  exp[-.693^)l‘-Jl)2]  (2-40) 

The  right  hand  column  vector  B  for  this  computation  was  the  self 
covariance  function  Cqq,  which  is  just  the  variance  of  the  electric 
field  on  the  near-field  measurement  plane.  (Similar  results  are 
obtained  if  a  measured  cross  covariance  function  is  used  as  the  known 
right  hand  column  vector).  The  computed  element  covariances  generally 
replicate  the  prominent  "lobes"  of  the  exact  element  covariances. 

Plots  of  the  magnitude  of  the  exact  and  computed  near-field  cross 
covariance  function  C^iq  (for  q”l,81)  are  presented  in  Figure  2-16. 
The  computed  curve  in  the  figure  agrees  closely  with  the  exact  curve 
over  the  portions  of  the  graph  where  C4iq  is  large.  The  close 
agreement  is  attributed  to  the  fact  that  functions  obtained  by  summing 
or  integrating  over  a  large  number  of  complex  quantities  are  usually 
insensitive  to  small  errors  in  the  summand  or  integrand. 

The  results  just  discussed  indicate  that  useful  estimates  of  the 
near-field  covariance  functions  can  be  obtained  from  a  knowledge  of  the 
self  covariance  function  on  the  measurement  plane.  Knowledge  of  at 
least  one  cross  covariance  function  is  also  needed  in  order  to  provide 
a  good  check  of  computational  accuracy.  Additional  research  is  needed 
in  order  to  determine  Che  feasibility  and  accuracy  of  this  method  of 
recovering  the  covariance  functions  for  arrays  comprised  of  hundreds  or 
even  thousands  of  elements. 

D.  Probability  Density  Functions 

A  thorough  characterization  of  the  radiation  pattern  statistics 
can  be  obtained  if  the  probal?ility  density  function  (p.d.f.)  for  either 
Che  radiated  field  amplitude  or  the  radiated  power  density  can  be 
discerned.  Once  the  p.d.f.  is  known,  the  standard  deviation  or  any 
other  higher  order  statistical  moment  can  be  computed.  Similarly,  the 
cumulative  probability  distribution  (c.p.d.)  is  readily  computed  by 
integrating  Che  p.d.f.  Of  course,  the  computations  of  Che  statistical 
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Figure  2-16.  The  exact  and  the  approximate  near-fleld  covariance 
function  for  the  electric  field  at  yO  with  the 
conjugate  of  the  electric  field  at  all  other  sample 
points  on  the  measurement  plane. 


tBoments  or  the  c.p.d.  mAy  have  to  be  performed  nuaericelly  because 
closed-fora  analytical  expressions  consisting  of  a  finite  sum  of 
elemetary  functions  for  Che  integrals  involving  the  p.d.f,  may  not  be 
available.  Nevertheless,  it  is  usually  more  efficient  to  compute 
radiation  pattern  statistical  data  via  the  use  of  a  reasonably  accurate 
p.d.f.  than  to  compute  radiation  pattern  statistical  data  directly  in 
terms  of  the  near-field  statistical  data  as  was  done  for  Che  average 
radiated  far  field  power  density. 

A  study  of  Che  modern  literature  concerning  electronagneCic 
scattering  and  radiation  statistics  indicates  that  the  p.d.f.  derived 
by  H.  Nakagami  [16]  in  1954,  and  chat  also  was  independently  de^ivf^d  by 
P.  Beckman  [17)  in  1961,  is  the  most  general  one.  Ocher  p.d.f. 's  for 
radiation/scattering,  such  as  these  due  to  S.O.  Bice,  Hoyt,  and  Lord 
Rayleigh  [  18-20  ]  auy  be  derived  from  the  Nakagami  p.d.f. 

The  key  assumption  in  Che  derivation  of  the  Nakagasu.  p.d.f.  is 
that  the  real  and  iswginary  parts  of  the  randomly  varying 
electromagnetic  field  both  follow  a  Gaussian  p.d.f..  The  Gaussian 
p.d.f.  is  an  excellent  approximation  to  the  true  p.d.f.  if  Che  electric 
field  is  produced  by  a  large  number  of  statistically  independent 
sources.  In  practice,  about  7  or  sure  statistically  independent 
sources  are  sufficient  to  yield  closely  Gaussian  p.d.f. 'a  for  the  real 
and  imaginary  parts  of  the  radiated  field. 

It  may  appear  at  first  sight  that  Che  Nakagami  p.d.f.  is  not 
applicable  to  the  randomly  excited  wire  array  studied  in  subsection  C 
because  the  element  currents,  as  well  as  the  measured  near-field 
electric  fields,  are  definitely  correlated.  However,  the  radiated 
electric  fields  at  any  point  in  the  near-field  or  the  far-field  of  the 
antenna  can  be  expressed  as  a  sum  of  terms  involving  non  random  matrix 
elements  multiplied  by  Che  statistically  independent  randoa  input 
voltage?.  Specifically,  the  electric  field  Eq  is  related  to  Che  input 
voltages  Vg  as 


E-  *  Z  2  G  r  V 
9  1  k  tk  tk'^k 


(2-41) 


where  all  other  symbols  are  as  previously  defined.  The  real  part  Xq  of 
Cq  and  the  imaginary  part  Tq  of  Eq  are  then 
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Xq  -  z  [{Hr)qp(Vr>p  -  {Hi)qp  (Vi)p  ]  and  (2-42) 

P 

Yq  -  £  [  (Hi)qp{Vr>p  ♦  (Br)qp(Vi)p  ]  (2-43) 

P 

wiiere 

(Vf)p  >  Che  real  part  of  Che  coaplex  random 
input  voltage  Vp, 

(Vi)p  “  the  imaginary  port  of  the  complex 
random  input  voltage  Vp,  and 

where  the  coefficients  of  (Vf.)p  and  (Vi)p  are 

(Hr)qp  -  Re  ,  and 

(Hi)qp  -  Im  {2,  } 

The  real  part  of  Eq  given  by  Equation  (2-42)  and  the  imaginary 
part  of  Eq  given  by  Equation  (2-43)  are  each  just  Che  sum  of 
staatiscical  ly  independent  random  variables,  and  if  the  nua^r  of 
voltages  is  i  7,  the  p.d.f.  for  X  and  the  p.d.f.  for  Y  will  both  be 
Gaussian  {  17  ]  to  good  approximation.  Hence,  Che  amplitude  of  Che 
radiated  field  will  follow  the  Nafcagami  p.d.f.  This  is  true  for  both 
Che  (radiating)  near-field  amplitude  and  the  far-field  amplitude. 

The  derivation  of  the  Nakagami  p.d.f.  involves  the  use  of  the 
joint  probability  density  of  X  and  Y.  The  joint  p.d.f.  for  two 
Gaussian  random  variables  X  and  Y  is  known  to  be  [  8  ]. 


(2-44) 

(2-45) 
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where  X  and  Y  are  the  average  values  of  X  and  Y,  and  aj  and  02 
are  the  standard  deviations  of  X  and  Y,  respectively,  p  is  the 
correlation  coefficient  for  X  and  Y,  and  v  ,  We  have  omitted 
the  subscript  q  frOM  X  and  Y  for  convenience. 

The  Nakagatni  p.d.f.  f(?>  for  the  power  density  P  is  preferred  here 
and  is  obtained  as  follows.  first,  the  joint  p.d.f.  for  X  and  Y  is 
expressed  in  polar  coordinates  by  replacing  X  and  Y  by 


X  =  Vp  cos(iJ>)  and  (2-47) 

Y  =  VT  sin(4i)  ,  (2-48) 

where  P  ■-  X^+Y^  ,  and  integrated  over  the  annular  ring  of  thickness 


lo  2-^  ,  to 

wit: 

dP  ] 

[  (if 

4iTaiO;^Y  j 

11  I  (v^  cos(4i)-<X>)^  .  f/'Psii.((^)-<Y>)^  I 


exp  I  2p{VF cps(|)-<X>HVF  sin(4.)-  <Y>]  j| 


(2-49) 

Next,  the  integral  is  expressed  as  an  infinite  susi  involving  Bodified 
Bessel  functions  of  the  first  kind  and  integer  order.  The  resulting 
expression  for  f(P)  is  then 


f(P) 


1 


exp  igP  +  hA  ] 


(?-;)]  •■'^1 


(2-50) 
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where  Cn  is 

Neuaanrf  s 

factor  (  En*l  for  n  “  0  and  en*^  otherwise)  and 

where  the  other  syabols 

are  defined  as  follows: 

&  ■ 

1  ,1  1, 

2  * 

(2-51) 

a  = 

Oj  +  02  + 

2  2 

[  4P^  +  (Oj  -  02)  ]  *S  , 

(2-52) 

e  “ 

Oj  +  Oj  - 

[  4p2  *  (oj  -  02)] 

(2-53; 

A  - 

<  X  + 

<  Y>^  , 

(2-54) 

a  ■ 

~  coa(«5i  -  dj)  , 

(2-55) 

b  » 

A  ain(52  -  ^2^  > 

& 

(2-56) 

^1  - 

-  tan''^ 

<X>  ] 

<Y>  1  ’ 

(2-57) 

1  tan  ^ 

2 

Oi-O^j  ' 

(2-58) 

h  ■ 

a  ^  cx  ^  b  ^ 

S  ,  and 

(2-59) 

c  - 

“"-‘[I 

[taa(«i-62) ] j  . 

(2-60) 

All  of  Che  new  quanCitiet  for  this  expression  for  f{p)  are 
conputed  solely  in  terns  of  the  sane  variables  that  are  needed  to 
define  the  joint  p.d.f«  for  X  and  Y.  These  are  (1)  the  average  values 
of  X  and  Y,  (2)  the  standard  deviations  of  X  and  Y,  and  (3)  the 
correlation  coefficient  for  X  and  Y,  It  should  be  noted  that  these 
variables  generally  depend  on  frequency,  array  scan  angle  and  the 
coordinates  of  the  field  point  under  consideration.  In  particular,  the 
power  density  P  ■  P(w,^  )  depends  on  frequency,  the  aziauth 
angle  of  the  field  point,  and  the  array  scan  angle. 

The  average  values  and  standaro  deviations  of  X  and  Y  as  veil  as 
correlation  coefficient  O  can  be  cosqiuted  with  the  aid  of 
Equations  (2-42)  and  (2-43)  when  the  statistics  of  the  input  voltages 
are  specified  froa  theory  or  expcriaecC.  The  equations  for  X  , 
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^  f  ^  It  O''  cooputed  directly  in  terms  of  the  input 

voltage  statistical  parameters. 

The  average  values  and  standard  deviations  of  X  and  Y  are  also 
readily  obtained  from  a  Type  II  near-field  measureaent  as  described  in 
subsection  C.  The  co(q>utation  of  the  correlation  coefficient  p  frosi 
the  measured  data  for  arbitrary  random  excitations  involves  problems 
similar  to  those  encountered  for  the  covariance  functions  that  vere 
discussed  in  the  preceeding  subsection.  However,  it  is  possible  to 
obtain  valid  estimates  of  P  by  assuming  that  the  random  voltage 
variations  have  the  same  variance  at  each  element. 

Mumerical  analyses  were  conducted  to  assess  the  validity  of  the 
Nakagami  p.d.f.  for  computing  the  far-field  amplitude  p.d.f.  for  the 
nine  element  wire  array  described  previously  in  subsection  C.  The 
computations  of  the  Nakagami  p.d.f.  were  accom;  lished  by  numerically 
integrating  the  expression  for  f(P)  given  by  Equation  (2-49)  rather 

than  by  suoaing  the  series  expression  given  by  Equation  (2-50).  These 

computations  were  performed  for  two  different  sets  of  input  parameters 
X  ,  Y  ,  oj ,  02i  and  p  to  the  Equation  (2-49).  The  first  set  of 

input  parameter  values  were  the  values  of  X  ,  Y  ,  Cj ,  02* 

computed  by  specifying  the  input  voltage  statistics  for  statistically 
independent  input  voltages.  This  set  forms  a  proper  set  of  input 
parameters  for  the  Nakagsmi  p.d.f.  The  second  set  of  values  for  X  , 
Y  t  °lt  ^2t  were  computed  directly  from  simulated  correlated 

voltages  with  an  assumed  Gaussian  correlation  function.  This  set  of 
input  parameters  involves  correlated  "sources"  and  therefore  violates 
the  strict  criteria  for  obtaining  the  Nakagami  p.d.f.  The  input 
voltages  for  both  situations  had  Gaussian  phase  distributions  with  zero 
mean  value  and  standard  deviation  of  57  degrees  on  all  elements.  The 
Gaussian  correlation  for  the  correlated  voltages  was 

-  exp  [ -[,693  ]  (2-51) 

where  ic  and  k*  are  element  numbers. 

Monte  Carlo  calculations  were  performed  to  produce  data  for 
constructing  the  cumulative  probability  distribution  for  the  far-field 
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amplitude  in  the  boresighc  direction  along  the  nose  of  the  aainbeam  and 
in  a  direction  of  a  low  sidelobe.  The  cumulative  distributions  were 
obtianed  from  30  Honte  Carlo  runs.  The  cuantlative  distributions  based 
on  the  Nakagaai  p.d.f.  were  obtained  by  nusierically  integrating  f(P). 

Plots  of  the  c\aaulative  distributions  are  shown  in  Figures  2-17 
through  2-20.  Each  figure  is  a  plot  of  the  probability  that  the  field 
amplitude  is  "less  than  or  equal  to"  the  abscissa.  The  abscissa  in 
Figures  2-17  and  2-19  is  the  relative  amplitude  in  real  numbers  while 
the  abscissa  in  Figures  2-18  «nd  2-20  is  the  relative  amplitude  in 
decibels.  The  sbscissas  were  scaled  in  this  Mnner  to  display  the 
interesting  fact  that  the  cumulative  distribution  resemble  Gaussian 
cuaulative  distributions  when  plotted  this  way,  in  agreement  with 
theory  for  the  mainbeam  and  low  sidelobe  regions  [  2  ].  Each  figure 
contains  two  cumulative  distribution  c’.rves.  These  are  (1)  the  curve 
constructed  from  Monte  Carlo  data,  (2'^  the  curve  obtained  from  f(P}. 
The  distributions  in  Figures  2-17  and  2-18  correspond  to  the 
statistically  independent  voltage  sources,  and  distributions  in  Figures 
2-19  and  2-20  correspond  to  the  correlated  input  voltages. 

Inspection  of  Che  plots  contained  in  Figures  2-17  through  2-20 
shows  that  the  Nakagami  curve  obtained  for  statistically  independent 
sources  is  in  very  good  agreement  with  the  Monte  Carlo  curve.  Exact 
agreement  is  not  expected  because  the  Honte  Carlo  curve  is  constructed 
from  "only"  50  runs.  The  Nakagami  curve  obtained  for  correlated 
sources  is  slightly  less  accurate  but  it  is  a  reasonable  approximation 
to  the  Monte  Carlo  curve. 

These  results  serve  to  verify  Naktgam's  p.d.f.  for  statistically 
independent  sources  and  they  also  indicate  Chat  it  may  be  a  useful 
engineering  spproxination  for  certain  classes  of  moderately  correlated 
sources.  Namely,  it  may  be  applicable  for  near-field  measurements 
involving  randomly  excited  waveguide  arrays  for  situations  where  the 
envelope  of  the  near-field  covariance  function  Cq'q  decreases  at  least 
as  rapidly  as  (1/  q-q'  )  as  the  difference  between  q  and  q*  becoaies 
"large."  This  condition  appears  to  be  crus  for  the  nesr-field  data  for 
the  wire  arrays  considered  in  this  study.  However,  further  studies  are 
needed  to  establish  Che  range  of  applications  for  idiich  the  Nakagami 


53 


Figure  2-17.  Cumulative  probability  distribution  for  the  radiated  power 
density  of  •  aalnbeam  or  grating  lobe  for  a  9-eIefflent  wire 
array  with  uncorrelated  random  input  voltages. 


RELATIVE  POWER  (dB) 

Figure  2~18.  Cumulative  probability  diacributlon  for  the  radiated  power 
density  of  a  low  aidelobe  for  a  9-element  wire  array  with 
uncorrelatad  random  input  voltages* 


* 


Figure  2-19.  Cumulative  probability  distribution  for  the  radla:ed  power 
density  of  a  malnbeam  or  grating  lobe  for  a  9-elenent  wire 
array  with  correlated  random  Input  voltages. 


Figure  2-20.  Cumulative  probability  distribution  for  the  radiated  powe 
den£.lty  of  a  low  aldelobe  for  a  9-element  array  with 
correlated  random  input  voltages. 


p.d.f.  can  furnish  valid  engineering  estiautes  ot  radiation  statistics 
for  correlated  sources  in  general  and  fcr  out-of-band  waveguide  arrays 
in  particular. 

A  sispler  but  less  accurate  approximate  expression  for  the 

radiation  aaiplitude  can  be  obtained  that  is  nevertheless  potentially 
useful  for  smse  engineering  applications  where  less  accuracy  can  be 
tolerated.  The  expression  is  derived  by  first  ignoring  the  correlation 
between  the  real  and  imaginary  parts  of  the  radiated  field.  The 
effects  of  near-field  source  correlstion  are  then  incorporated  as  part 
of  a  paraaieter  that  appears  naturally  in  the  expression  that  was 

derived  by  ignoring  the  correlation  between  the  real  and  imaginary 
parts  of  the  radiated  field.  The  approximate  expression  for  the 

radiated  power  density  obtained  in  this  manner  is 


1 


P(t).4i)  4- 

x(tis) 


•  I  [B(ui.4>)3 

o 


(2-62) 


where  f[  P(w  ,(fi  )]  “  probability  density  function  for  the  far-field 

power  density, 


(2-63) 


tCoj)  =  r,  E  C  ,  exp[j-  sin(iji)(Y  ,-T  )1  , 
q  q  q  q  c  q  o 


(2-64) 


I  [  ■]  B  Modified  Bessel  function  of  the  first  kind 
and  order  zero,  and 


(2-65) 


(2-66) 


~  T(a;) 
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Plots  of  the  cusKilative  distribution  coaputed  frooi  this  equation 
are  compared  with  the  cuawlative  distribution  obtained  from  the  Monte 
Carlo  data  are  shown  in  Figures  2~21  and  2-’ 22.  Resonably  good 
agreement  is  obtained  between  the  Monte  Carlo  curves  and  the  curves 
computed  from  Equation  (2-62).  However,  additional  testing  is  needed 
to  properly  define  the  limits  of  applicability  of  this  relatively  crude 
p.d.f.  for  antenna  problems. 

E.  Tisw  Domain  Statistics 

It  is  a  well  known  fact  that  the  time  domain  response  of  a 
deterministic  radiating  system  can  be  analysed  conveniently  from  a 
knowledge  of  the  complex  frequency  response  via  the  Fourier  transform 
technique.  The  same  procedure  can  also  be  employed  to  obtain  the 
statistical  average  tine  domain  response  of  a  randomly  excited 
radiating  system.  The  probability  density  function  (p.d.f.)  for  the 
time  domain  response  can  also  be  discerned  from  the  stochastic 
equations  obtained  in  this  manner. 

The  behavior  of  the  far-field  electric  field  as  a  function  of  time 
t  is  obtained  from  Fourier  analysis  of  the  frequency  response  as 

E(t,.*  )  *  :^  I  E(w  .♦)exp(iu  t)  (2-67) 

n  n  n 

where 

E(t,i^  )  >  complex  electric  field  at  azimuth  angle  ^  versus  time, 
and 

E(<»ijj,t)  *  the  complex  electric  field  at  azimuth  angle  ♦  for 

radian  frequency  u  "Zst, where  f  is  the  frequency  in 
Hertz. 

E(  <1 )  is  the  product  of  the  input  pulse  spectrum  H(  io^)  and  the 
frequency  response  E®(  ♦ )  obtained  from  a  uniformly  weighted  input 

frequency  spectrum  of  width  exceeding  at  least  twice  the  bandwidth  of 
Hfi^ij).  Hence, 

E(  )  -  E«»(  •  H(  (j^)  (2-68) 
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Figure  2-21.  Cumulative  probability  distribution  for  the  radiated  power  density 
of  the  malnbeain  or  grating  lobe  for  a  9-eLeinent  array  with  un¬ 
correlated  input  voltages. 


L 


The  power  density  P(t,$  )  is  then  coaputed  as  E*(t,  ♦)  •  E(c,  and  i# 
written  explicitly  as  the  double  sub 


P(t,  ^ > 


^n  “n"  *  “n*  *^****  ^  ^^  “’n”  ^  (2-68) 


The  tine  doaain  behavior  of  a  pulsed  stochastic  antenna  systea  can 
be  obtained  with  the  aid  of  Che  preceding  f requency-doBain  equations 
and  the  Fourier  transfora  with  respect  to  frequency.  In  particular, 
the  tine  dependent  atacistical  average  conplex  far-field  electric  field 
may  be  written  as 


<E(C,<>»  “  ^  exp[j(o^t] 


(2-69) 


In  analogy  with  Che  techniques  described  previously  for  coBpuring 
average  power  density,  the  statistical  average  time  dependent  power 
density  is  written  as 

<1*(t,  <E* (!,♦)>•  <E(t,<t.)>  (2-70) 

^n  ^n  ^n’^n  ^n'^a  < V  ' 

(n’fl  n) 

where  C  standard  deviation  of  E(  is  ,^  ), 
n  n 

«  standard  deviation  of  E*((i>  )»  snd 

■  cross-correlation  function  of  E(u>^,^  )  and  E*(ft)^,,i>  ), 

The  cross-correlation  function  ^qI^(  *^ )  appearing  here  is  not  a  tijie 
correlation  function  but  rather  is  a  freqjcncy-depeudent  correlation 
function  for  the  electric  field  and  its  conjuga^.v^  al  different 
unless  the  radiating  systea  contains  non-linear  devices  that  have 
pronounced  hysterisis  or  contains  devices  whose  characteristic  response 
time  is  coBoensurate  with  the  input  pulse  width. 
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The  p.d.f.  f{  P(t,  ^)}o£  Che  tine  doMin  power  density  P(t,  #  )  is 
the  Nakagani  p.d.f.  when  the  field  and  its  conjugate  at  different 
frequencies  are  uncorrelated.  Thip  can  be  deduced  at  once  by 
recognising  that  the  reel  and  imaginary  parts  of  the  electric  field 
E(t,4  )  each  consist  only  of  suns  of  statistically  independent  terns 
when  the  Are  sero.  The  input  paraneters  for  the  Nakagani  p.d.f. 

are  the  average  values  of  the  real  and  isuiginary  parts  of  }, 

their  standard  deviations,  and  the  correlation  coefficient  for  the  real 
and  imaginary  parts  of  E(t,  ^).  The  correlation  coefficient  for  the 
real  and  isuginary  parts  of  E(t,  #)  is,  of  course,  non-zero  even 
when  the  frequency  cross  correlation  coefficients  R  ,  are  zero. 

A  nunterical  aiuilaticn  of  the  pulse  distortion  suffered  by  a 
rectangular  pulse  radiated  by  a  reflector  antenna  for  both  in-band  and 
out-of-band  frequencies  for  deterministic  conditions  is  presented  and 
discussed  in  Section  V. 
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SECTION  Ill 

NEA&-FIELD  ANTENNA-ANTENNA  COUPLING 


A,  Introduction 

An  exact  analysis  of  Che  coupling  between  real-world  cosited 
antennas  would  be  extrenely  difficult  or  perhaps  iapossible  to  achieve. 
Accordingly,  the  efforts  on  this  task  were  devoted  to  deriving 
approxisate  but  accurate  coupling  prediction  equations.  Three 

different  techniques  for  coupling  analysis  were  considered.  These 
three  techniques  are  denoted  respectively  as  (1)  Plane  Wave  Spectrun 
(PWS),  (2)  Spherical  Spectrusi  Wave  (SWS),  and  (S)  the  Ceoaetrical 
Theory  of  Diffraction  (CTO).  Most  of  Che  research  efforts  were  devoted 
toward  development  of  the  PWS  technique.  Multiple  scattering  effects 
are  not  addressed  in  Che  analysis.  However,  aultiple  scattering 
effects  can  be  approximately  analysed  if  the  scattering  siatrix  of  each 
antenna  is  specified  from  theory  or  measur«Bents .  Theory  and  equations 
for  the  PWS  and  SWS  approaches  are  presented  and  discussed  in 
Subsection  B  and  C,  respectively.  Soim  results  of  numerical 
simulations  for  PWS  approach  are  also  presented  in  Subsection  B.  The 
GTD  technique  for  deterministic  antenna  analysis  has  been  adapted  to 
yield  comparatively  aisqtle  equations  for  certain  classes  of  antenna 
coupling  problems.  Discussions  of  the  GTD  technique  are  contained  in 
subsection  D. 

The  analyses  presented  in  this  section  were  conducted  in  the 
frequency  domain.  The  primary  goal  was  to  obtain  equations  vslid  over 
wide  frequency  bandwidths.  The  tisK  doauin  response  can  then  be 
obtained  by  nusKrically  computing  the  Fourier  Transform  of  the 
frequency  domain  equations.  The  wideband  frequency  response  of  swept 
CW  radiating  systems  can  be  characterised  by  numerically  computing  the 
frequency-averaged  pattern  and  standard  deviation,  as  described  in  the 
previous  section. 

B  .  PWS  Analysis 

The  general  plane  wave  apectruo  scattering  matrix  theory  was 
originally  fcnsulated  to  treat  antenna  coupling  problas  and  it  is  well 
documented  in  Reference  [  21} .  The  theory  and  equations  sre  directly 


applicable  to  antenna  coupling  situations  shown  in  Figure  3-1  where  one 
antenna  is  in  the  forward  half-plane  of  the  other  antenna  and  the  relative 
rotation  angles  a  and  B  defined  in  Figure  3-2  are  leas  than  about  70 
degrees  in  aagnitude.  The  application  of  Che  FWS  analysis  Co  other 
arrangeaenCs  that  do  not  obey  these  criteria  depends  on  the  ability  to 
select  suitably  "canted"  reference  planes  with  appropriate  spatial  filters. 
This  is  a  recoomended  area  for  future  research.  The  FWS  analysis  for  the 
canonical  situation  depicted  in  Figure  3-1  is  discussed  in  the  following 
paragraphs . 

The  complex  voltage  V(R,  induced  in  Antenna  B  when  illuminated 

by  Antenna  A  is  derived  from  FWS  analysis  as 


p>  ~  ^pi^^zn  ~  ^*q)*exp[  “jlt«n  * 

(3-1) 

%rtiere 

A(ky^,k2n)  "  transverse  vectorial  plane  wave  spectra  of  Antenna  A, 

B(kyB  ~  ^yp>  ~  *  shifted  transverse  vectorial  plane 

wave  spectra  of  Antenna  B, 

hy*  *  mt*»  value  of  ky  in  the  range  ■*  ko  *'^y  - 

k^n  •  n^**  value  of  k^  in  the  range  -  kg  ^  k^  k^, 

syp  “  kosinCctp}!  evaluated  at  the  specified  azisutb 
rotation  angle  Op, 

^zq  ”  Vo*^“(®qK  evaluated  at  the  specified  elevation 
rotation  angle  Bq, 

*  ^xmn  *  *^ym  7  *»  where  kMn  “  (kj  -  ky^  ~ 

and  Che  x,  y  and  z  are  unit  vectors, 
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R  =Xx  +  Yy  +  ZZj  where  X,  Y,  and  Z  are  the  cartesian 
coordinates  of  the  center  of  rotation  of  Antenna  B 
relative  to  the  center  of  rotation  of  Antenna  A,  and 

Cq  =  a  frequency-dependent  proportionality  factor. 

The  proportionality  factor  Cq  is  independent  of  antenna  arrangements  and 
orientations  and  will  be  suppressed  in  the  .-emaining  equations.  The  power 
coupled  in  Antenna  B  is 

otp,  3q)  =  !v(R,  “p,^q)  'UO.S)  *  Tj.]  (3-2) 

p  , 

where  ‘  j.  is  the  real  part  of  the  complex  admittance  at  the  output  port. 

The  PWS  of  Antenna  A  may  be  expressed  in  terms  of  the  propagating 
modes  in  the  waveguide  feed  as 

ACkyjjjk^n)  ”  ^ 

(3-3) 

where 

a^  =»  the  complex  excitation  coefficient  of  the  ^th 

waveguide  propagating  mode, 

^(Cj^^)  =  optical  transform,  computed  as  (1/2)  •[!  +  cos 

H^(Cjj^)  ••  the  far-field  vectorial  pattern  of  the  feed  horn 
when  only  the  th  mode  propagates , 

=  points  on  the  reflector  aperture. 

The  angle  and  the  reflector  apertire  coordinates  are  depicted  in  Figure 
3-3.  Equation  (3-3)  can  be  rewritten  as  the  sum  of  the  reflector 
transverse  vectorial  spectra  A  ^  (kyg,,k2i,)  , 

A(kyjj,k2n^  ■  Z^a^  A^(kyni,k2n)  •  (3-A) 

The  corresponding  equations  for  the  PWS  of  Antenna  B  are 
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(3-5) 


the  alterr.::tive  form, 


B(kynj,kj5n)  “  J,  b  B  (k  ,k  ) 
*  PUP  ym  zn 


(3-6) 


where 


b  =  the  complex  excitation  coefficient  of  the  waveguide 
propagating  mode,  and 


all  other  symbols  are  as  previously  defined. 

The  equation  for  the  statistical  average  power  P(R,ap,Bq) 
coupled  to  Antenna  B  when  illuminated  by  Antenna  A  is  derived  as 

\  {K’>  <%>  *  'V 


—A  — r  ^  1 

X  A  ,  (k  ,k  ,)*B  ,(k  ,  -  k  ,k  ,  -  k  )exp  jk  ,  , 'R 

p'  ym’  zn'  p’  ym'  yp  zn '  zq  ^  m'n'  J 


X  A  (k  ,k  )*B  (k  -  k  ,k  -  k  )exp 
P  ym  zn  p  yn  yp  zn  zq 


where  the  statistical  quantities  are  defined  as 


H'«-  'l| 


(3-7) 


<^a  >  ■  average  value  of  the  complex  mode  excitation 

^  coefficient  a^  , 

-  average  value  of  the  complex  mode  excitation 
coefficient  , 

C^,  “  covariance  function  for  a*  and  a  .  and 

P  P  p  p 

C^i,  =  covariance  function  of  b*  and  b 

P  p  p  p 

The  wave  vector  k„n  is  defined  as  kn,n  -  k^n  ^^yn  Y  * 
km'n'  is  defined  as  lqn«R'  *  kxm'n'  *  +  ^'ym'  Y  +  ^zn 


Similarly, 
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Nunerical  sinulations  were  conducted  to  coapute  the  power  coupled 
between  a  pair  of  cosited  antennas.  The  antenna  arrangesKnt  is  shown  in 
Figure  3-4.  Antenna  B  is  displaced  a  longitudinal  distance  X  and  Lateral 
distance  Y  frost  Antenna  A.  The  antennas  are  not  displaced  in  the  vertical 
Z  direction.  The  power  coupled  to  Antenna  B  versus  the  rotation  angle  was 
computed  for  selected  values  of  x  and  y  for  rotation  angles  in  the  range  - 
30  degrees  to  >30  degrees.  Only  the  azimuth  Plane  Wave  Spectrums  of  the 
antennas  were  used  in  the  calculations. 

The  apertures  of  Antenna  A  and  Antenna  B  are  24  inches  and  48  inches, 
respectively.  The  antenas  are  fed  by  WR-187  waveguide  whose  cutoff 
frequency  for  the  TE20  mode  is  6.309  GHz.  The  in-band  operating  frequency 
of  each  antenna  is  5.5  GHz.  Calculations  were  made  for  the  in-band 
frequency  of  5.5  GHz  and  the  out-of-band  frequency  of  6.5  CRs.  The  power 
flow  in  each  waveguide  feed  is  entirely  in  the  TEjq  mode  at  5.5  GHz.  The 
TE|o  and  TE20  modes  can  both  propagate  at  6.5  GHz. 

Antenna  coupling  can  be  conveniently  described  in  terms  of  the  mutual 
gain  M(R,a)  relative  to  a  pair  of  isotropic  radiators.  M(R,  a)  is  computed 
as 

M<R,  a)  -  ,  (3-8) 

'  '  a 

where  ^  is  the  operating  wavelength  and  Fj,  is  the  power  input  to  Antenna 
A.  For  the  purpose  of  this  study,  a  power  transfer  efficiency  factor  of 
1.0  is  assumed  for  Antenna  A  and,  hence,  the  total  power  radiated  by 
Antenna  A  is  assumed  to  equal  the  input  power. 

The  in-band,  far-field  antenna  patterns  of  Antenna  A  and  Antenna  B  are 
plotted  in  Figures  3-5  and  3-6,  respectively.  The  estimated  theoretical 
gain  relative  to  an  isotropic  radiator  Is  indicated  at  the  top  of  each 
plot.  The  antennas  were  regarded  as  linear  apertures  with  an  "aperture 
efficiency''  of  0.9  relative  to  comparable  uniformly  illtaainated  apertures 
for  purposes  of  estimating  the  gain. 

Out-of-band  antenna  patterns  for  the  two  antennas  are  shown  in  Figures 
3-/  and  3-6.  Figure  3-7  shows  the  o«it-of-band  pattern  for  Antenna  A  for 
equal  power  flow  in  both  swdes  with  the  phase  of  the  TE20  mode  >35° 
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Sketch  of  the  near-fleld  arrangement  of  Antenna  A  and 
Antenna  B  used  In  the  numerical  simulations. 


dB/lSO 


of  5.5  GHz  for  wav«guide  feed  power  flow  entirely  in 


16.1  dB/lSO 
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frequency  of  6.5  GHf.  for  equal  waveguide  power  flow  in  the 
TE^-  and  TE-  modes  with  relative  phase  angle  of  35  degrees. 


relative  to  the  TEjq  mode.  Similarly,  Figure  3-8  shows  the  out-of-band 
pattern  for  Antenna  B  for  equal  power  flow  in  both  modes  with  the  phase  of 
the  TE20  mode  -50°  relative  to  the  TEjq  mode. 

The  in-band  to  in-band  mutual  gain  of  the  antenna  pair  is  shown  in 
Figures  3-9  and  3-10.  Each  figure  shows  a  plot  of  the  mutual  gain 
expressed  in  decibels  versus  the  rotation  angle  o  in  degrees  for  the 
indicated  values  of  longitudinal  separation  distance  X  and  transverse 
separation  distance  Y.  The  mutual  gain  relative  to  a  pair  of  isotropic 
antennas  is  displayed  at  the  top  of  each  plot. 

The  out-of-band  to  out-of-band  mutual  gain  of  the  antenna  pair  is 
shown  in  Figure  3-11  for  the  indicated  out-of-phase  modal  propagation  for 
the  longitudinal  separation  distance  of  20  feet  and  lateral  displacement  of 
5  feet. 

Two  general  trends  common  to  the  in-band  to  in-band  and  the  out-of- 
band  to  out-of-band  situations  may  be  discerned  from  inspection  of  the 
plots  shown  in  Figures  3-9  through  3-11  and  from  additional  plots  contained 
in  Reference  [5].  First,  the  peak  value  of  the  mutual  gain  decreases  with 
increasing  lateral  displacements  for  a  fixed  longitudinal  displacement. 
Second,  tha  peak  value  of  the  mutual  gain  increases  with  increasing 
longitudinal  displacement  for  a  fixed  lateral  displacement.  Both  of  these 
trends  are  consistent  with  theory.  However,  the  peak  mutual  gain  may 
exhibit  "peaks"  and  "valleys"  if  the  computation  were  made  for  a  "finely- 
grained"  set  of  coordinates.  Nevertheless,  the  "envelope"  of  the  peak 
mutual  gain  would  exhibit  the  cited  trends. 

The  peak  mutual  gain  for  the  in-band  to  in-band  situations  always 
occurs  at,  or  very  near  to,  the  rotation  angle  at  which  Antenna  B  points  at 
the  center  of  Antenna  A,  whereas  the  peak  mutual  gain  for  the  out-of-band 
to  out-of-band  situations  is  shifted  a  few  degrees.  The  direction  and 
magnitude  of  the  angular  shift  depends  on  the  mode  excitation  in  the  feeds 
of  of  the  two  antennas.  The  magnitude  of  the  peak  mutual  gain  also  varies 
with  mode  excitations.  All  of  these  trends  are  consistent  with  theory  in 
the  sense  of  "coarse-grained"  behavior  mentioned  previously. 

Next,  consider  the  situation  depicted  in  Figure  3-12  where  the 
transmitting  antenna  illuminates  a  cylindrical  obstacle  located  between  the 
transiultting  and  receiving  antennas.  In  order  to  determine  the  effect  of 
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4  dB/lSO 


(0p)  NIV9  ivninw 


distance  Y  ■  5  feet 


the  interposed  obstacle  on  the  antenna  coupling  the  total  transverse 
spectrua  for  the  transmitting  antenna  and  Che  obstacle  which  is  incident 
upon  Che  receiving  antenna  saist  be  detendned.  This  calculation,  which  is 
fully  described  in  References  [21]  and  [22],  is  Che  basis  for  the 
evaluation  of  the  far-field  pattern  degradation  caused  by  the  near-field 
obstacle.  For  the  case  of  Che  interposed  obstacle,  Che  spectrua  incident 
on  Che  receiver  is  Che  sua  of  the  CrsnsaitCing  antenna  spectrua  A(k)  and 
the  obstacle  scattered  spectrua  A*(k)  given  by  Equation  (5-8)  of  Section  V. 
Assuaing  that  the  aulciple  reflection  effects  are  negligible,  the  ccmplex 
voltage  response  can  be  coaputed  froa  Equation  (3-1)  by  replacing  the 
antenna  spectrua  A(k)  with  Che  coaposite  spectrua 

At(k)  «  A(k)  +  A  (k)  (3-9) 

8 

The  aethod  of  coaputstion  of  the  integral  is  thus  Che  saae  whether  or 
not  an  obstacle  is  present.  The  only  difference  is  whether  the  antenna 
spectrua  A(k)  or  the  sua  of  Che  antenna  and  obstacle  scattered  spectrua 
At(k)  is  employed  for  the  coaputation. 

The  plane  wave  spectrua  scattering  analysis  algorithm  was  used  to 
construct  a  data  file  corresponding  to  the  sua  of  Che  transmitting  antenna 
and  obstacle  scattered  spectra.  A  second  data  file  is  then  loaded  with  the 
receiving  antenna  spectrua,  and  the  transmission  integral  evaluation  is 
performed. 

Recall  that,  for  a  call  cylindrical  obstacle,  the  scattering  matrix 
eleaenC  function  in  the  vertical  (^  ■  constant)  plane  is  of  Che  fora  SIN 
(Z)  /(Z).  The  antenna  spectra  in  Che  vertical  direction  were  also 
approximated  by  appropriately  chosen  functions  of  SIN  (Z)  /(Z)  fora.  These 
functions  were  chosen  to  produce  Che  correct  3-dB  beaawidth  for  the  test 
antennas.  The  coaputation  can  be  further  simplified  if  the  coupling 
integral  is  performed  only  over  the  aziaSiCh  plane  (i.e.,  6  *  x/2,  ky 

variable).  This  approximation  is  termed  a  "linear  spectrua"  approxiaation, 
and  it  was  found  Chat  this  approximation,  for  the  cases  considered  herein, 
was  a  good  approximation  to  the  use  of  the  complete  (kyik^)  spectra.  Thus 
the  "linear  spectrua"  approxiaation  was  used  to  derive  the  antenna 
decoupling  data. 
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Measured  data  foe  both  6-inch  diaaeCar  and  24-inch  diasecer  cylinders 
located  between  48-inch  diaaeter  paraboloidal,  vertically-polarized 
reflector  antennas  operating  at  5.5  GBa  were  available  froa  previous 
aeatureaent  pregraaa  [  23  ]  .  Figure  3-12  shows  the  antenna/obstacle 

geoaetry.  The  calculated  antenna  decoupling  caused  by  the  interposed  6- 
inch  and  24-ioch  diaaeter  cylinders  for  antenna  separation  distances  of  17- 
feet  and  47~feet  is  shown  in  Figure  3—13.  The  term  "decoupling"  is  defined 
as  the  ratio  cf  the  received  voltage  with  the  obstacle  present  to  the 

received  voltage  with  the  obstacle  absent.  Thus,  the  decoupling  represents 

the  decrease  in  mtual  coupling  due  to  the  interposed  obstacle. 

As  shown  in  Figure  3-13,  the  calculated  and  aest  ured  results  for  both 
the  6-inch  and  24-inch  diaaeter  cylinders  for  the  47-foot  separation 

distance  are  in  very  good  agreeaent  with  the  average  aeasured  data.  We 

note  that  Che  actual  Measured  data  exhibit  a  periodic  oscillation  about  the 
average  values  presented  herein  due  to  aultiple  reflections  between  the 
obstacle  and  antenna.  Since  these  reflections  are  not  considered  in  the 
present  analysis,  it  is  appropriate  Co  coapare  the  calculations  with  these 
average  data.  The  typical  aaxiaua/ainiaua  bounds  on  the  aea&ured  data 
deviate  froa  the  average  values  by  approxiaately  *1.5  dB  and  *0.5  dB  for 
the  24-inch  diaaeter  and  6-inch  diaaeter  cylinders,  respectively. 

The  calculated  decoipling  for  the  17-foot  separation  distance  for  the 
6-inch  diaaeter  cylinder  is  also  in  very  good  agreeaent  with  the  Measured 
data.  However,  Che  results  for  Che  17-foot  separation  distance  for  the  24- 
inch  diaaeter  cylinder  exceed  the  Measured  data  by  a  few  decibels  due  to 
the  spproxiaationa  chat  were  aade  in  Che  analysis. 

Plots  of  the  Mssured  and  coaputed  nutual  gain  versus  the  azimuth 
pointing  angle  of  the  receiving  antenna  are  shown  in  Figures  3-14(a>  and  3- 
14(b}.  The  mutual  gain  plots  shown  in  Figure  3-14(b}  were  obtained  with  no 
obstacle  present.  The  Mutual  gain  plots  shown  in  Figure  14(b)  were 
obtained  obtained  with  a  12.(>-foot  Call,  2.0-foot  diaaeter  mast  located 
aidway  between  the  two  4.0-fc<ot  diaaeter  paraboloidal  reflector  antennas 
that  were  separated  17.0  feet  apart.  The  coaputed  data  were  obcained  vis 
Equations  3-1  and  3-8.  The  elevation  patterns  of  Che  antennas  were 
characterized  in  the  approximate  manner  described  previously,  and  multiple 
reflections  were  ignored.  Nevertheless,  good  agreement  between  measured  and 
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12.0-foot  tall,  4. 0-foot  diameter  mast  located  mldviay  (X  “8.5  feet)  between 
the  two  antennas  for  the  indicated  antenna  separation  diicance  R*17.0  feet, 
vertically  polarized  signals  at  ,'.5  GHz.  The  aperture  diametar  of  each 
paraboloidal  reflector  antenna  is  4.0  feet. 


coApuced  data  is  achieved  for  this  saaple  situation  as  well  as  for 
situations  involving  cylindrical  obstacles  of  O.S-foot  diavier.er  and 
4.0-foot  diaaeter. 

Application  ol  the  PWS  equation  to  predict  mstual  coupling  with  an 
obstacle  located  between  the  two  antennas  shown  in  Figure  3-1  and 
displaced  both  vertically  and  laterally  off  the  line  joining  the 
antenna  rotation  centers  is  well  understood.  The  anttnna/obstacle 
situation  just  . jscribed  can  be  analyzed  via  existing  equations. 
However,  the  application  of  PWS  analysis  to  handle  the  effects  of  an 
intervening  obstacle  located  between  the  two  antennas  depicted  in 
Figure  3-15  or  a  non-intervening  obstacle  such  as  a  ccsioion  iaouati;ig 
pole  to  which  two  or  nore  antennas  are  attached  are  areas  where  further 
research  is  needed. 

C.  Spherical  Wave  Spectrua  Coupling  Analysis 

The  formulation  of  the  SWS  coupling  analysis  can  be  illustrated 
with  the  ard  of  the  cosited  antenna  pair  ahosm  in  Figure  3-15  and 
the  antenna  coordinate  system  depicted  in  Figure  3-16.  The  various 
symbols  contained  in  Figures  3-lS  and  3-16  are  defined  as 

R  »  vector  from  center  of  Antenna  A  to  center  of  Antenna  B, 

R'  vector  from  the  center  of  Antenna  A  to  a  specified 
near-field 

R-R'  ^  vector  from  the  center  of  Antenna  t  to  the  specified 
near-field  point, 

a  •  azimuth  rotation  angle  between  Antenna  A  and  Antenna  B,  and 
~  elevation  rotation  angle  between  Antenna  A  and  Antenna  B. 

The  complex  voltage  response  V(R,  B  ,  a  )  of  Antenna  B  when 
illuminated  by  Antenna  A  can  be  written  as 

V{R,  3  ,  o)  -  ^  (R-R')  1®  ,  »  >  •  E*  (R'  >  dS'  (3-10) 

S* 

v:here  is  a  frequency  dependent  factor  and  where  and  denote  the 
complex  near-field  electric  fields  of  Antenna  A  and  Antenna  d. 


respectively.  The  surface  integral  in  Equation  (3-10)  is  integrated  over 
the  surface  of  the  sphere  of  radius  R’ . 

Tfie  vertically-polarized  and  horizontally  polarized  coBq>onen\’.s  of  the 
electric  field  of  Antenna  A  on  che  sphere  of  ra<‘ius  H'  centered  at  Anteina 
A  nsay  be  co-aputed  as 

tR',  0’,  4>')  “  I  a''  h^^^(KR')L  (cos  e')e"^“'^'  (3~11) 

V  ID  n  flm  A  nn 

for  vertica'  polarization  and 

Ef(r.*.  0',  <J.')  -  Z  Z  A^  h^^^{KP' )L  (cos  6’)e'-’*  (3-12) 

n  '  9  n  so  XX  tm 

(2) 

for  horizontal  polarization.  The  syi^ol  h^  denotes  the  spherical  Hankel 
function  'f  the  second  kind  and  order  n,  and  the  symbol  detotes  the 

associated  Legendre  polynomial  of  the  first  kind  of  order  n  and  degree  m. 
The  coefficients  and  A^  are  readily  computed  if  E®  is  known  on  any 

spherical  surface.  It  is  assumed  that  the  vertically-polarized  and 

horizontally-polarized  far-field  patterns  are  knotra,  and  consequently,  the 
coefficients  can  be  computed  from  them.  Equations  analogous  to  Equatiors 
(3-11)  and  (3-12)  can  he  written  for  the  vertically-polarized  and 

horizontally-polarized  fiel.*>  of  Antenna  E,  and  a  coordinate  transformation 
can  be  applied  to  refer  the  fields  of  Antenna  B  to  the  Antenna  coordinates 
or  vice,  versa. 

The  power  coupled  into  Antenna  B  due  to  illumination  by  Antenna  A  is 

then 

?(R,  d  ,  u)  Atm  *®n*m'  t o)  (3—13) 

n'm 

where  A  =  a'^  6  +/.**#  and  B  «  B^  0  B^  ^  .  The  coefficients  A  and  B 

and  the  integral  I(R, B  ,a}  are  computed  as 

ir 

F®  (  5,  <}•)  *^“®‘*®‘**  »  (3-U) 
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F^(6,*)L  ,  ,(cos8)e  sinOddd^,  and  (.1*15) 

A  n 


1  (R.P.a)  -  R 

na 


2ii  f 

vr  r  »'n{KR')hn(Kli-R' i>lTO»'‘^o»9’>I^’a’^co8(a'. 


e-j«i(^'-a:<s.;rt(2«  )de'd4,' 


(3-16) 


“  fl 

The  quantities  F  (0,<^)  and  F  (t,,  v)  are  the  total  vectorial  angular 
spe  truna  of  Antennas  A  and  B,  respectively.  These  total  spectruau  are 

""fll  ~ 

just  the  sua  of  the  vectorial  epsetruns  t  or  for  each  propagating  node 
of  the  antenna  feed  systoa.  Hence, 


F^'ve,^)  -  f“  (6,*)  and 


(3-17) 


(3-18) 


where  a^  and  denote  coaplex  excitation  coefficients  for  the  propagating 
waveguide  aodes. 

The  statistical  average  power  is  obtained  In  the  straightforward  but 
tedious  procedure  described  previously.  The  resulting  expression  is 


I  Z  Z  ,Z  ,Z  Z  I  ,l  ,  <A  •  A  Xb 
nmn  ot  pqp  q  ^nia  pq"^  ^  n  m*  p  q  ' 


I  (R,B,u)l  (g.B.a) 
nn  pq 


(3-19) 


The  joint  statistical  aocMot  ^A^  ’  is  coatputed 


9i 


2ti  u  2 It  it 

<^.‘:  >■//// 


<F‘^{9,<ti)-  L  (eose)L  {cos6'> 

ao  pq 


o  o  o  o 


(3-20) 

X  expf  j (p-n)<!)  1  8ln(6)  sin(0')  d0'  d^'Oddd^ 


and  •  H*p«q'!>  is  oomputed  ajjailarly.  Tbe  joint  statistical 

moiBent  «CF®(9  ,<i  )  •  (  P'  ,  »  is  obtained  from  a  knowledge  of  the 

waveguide  tuodal  coefficient  excitation  statistics,  to  wit: 


<Faie,((.)*p**(e' ,({,*)>  - 


I  z 

u  P 


,<a  a*>F^{e.*K^*{0’,(|)’) 
y  y  y  y  9^  ' 


(3-21) 


and  similarly  ior  <CF^(0  ,<}•  ) ‘F**  (8  ‘  ,^> '  )>  . 

Finally,  the  joint  statistical  ncments  <a^ 

P*  must  be  specified  from  theory  or  experiment.  the  deterministic 
analysis  is,  of  course,  recovered  from  the  statistical  equations  in  the 

limit  of  Vanishing  covariances  C  ,  for  all  indices  p  and  p*  . 

P  P 

The  rigorour  spherical  wave  spectr;ai  analysis  is  considerably  more 
involved  than  the  plane  wave  spectrum  analysis.  However,  the  spherical 
wave  spectrum  technique  has  the  very  desirable  attribute  of  being 
applicable  to  arbitrary  ar;.argemencs  and  orientations  of  a  near'field 
antenna  pair.  Nonetheless,  computer  run  times  may  be  excessive  for 
routine  use  of  the  SUS  technique  in  its  present  "rigorous"  form. 
However,  it  appears  that  a  substantial  reduction  in  the  computer  run 
time  may  be  achieved  without  seriously  degrading  the  accuracy  of  the 
SWS  analysis  by  derivins  approximate  jeries  representations  of  the 
integrals  via  cither  the  residue  calculus  or  replaceaent  of  the  full 
spherical  wave  series  by  asymptotic  series.  Moreover,  it  appears  that 


,  r 
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further  theoretical  and  ntaKrical  studies  of  the  SWS  technique  would 
ex;Midice  the  development  of  valid  WS  equations  for  analyzing  arbitrary 
antenna/obstacle/antenna  geometries.  Accordingly,  the  technique  is 
considered  to  be  a  potentially  valuable  area  for  candidate  further 
study  in  future  investigations. 

_ GTD  Analysis 

Relatively  simple  equations  are  presented  in  R.'fercnce  30  for 
estimating  the  coupling  of  power  between  cosiced  reflector  antennas. 
The  equations  were  derived  with  aid  of  the  GTD  analysis.  They  are 
applicable  to  situatinns  whereby  the  "mainbesas"  of  the  two  antennas 
ars  noL  directed  toward  each  other.  That  is,  the  equations  are 
applicabla  to  "sidelohe-co-sidelobe"  coupling  situations.  They  are  not 
anplicable  Co  "Kainbeas-to-maipbeasi''  or  ''mainbeam-to-sidelobe"  coupling 
situations.  Thus,  the  STD  equations  complement  the  current  WS  method 
presented  previously  which  is  particularly  well  suited  for  the 
flainbeam~tO~aainbeaa  and  mainbeam-to-sidelobe  coupling  situations.  It 
should  be  noted  that  the  equations  lack  the  rigor  associated  with  the 
PUS  and  SWS  analyses  and  that  it  is  difficult  to  judge  the  absolute 
accuracy  of  the  equations,  nevertheless,  the  equations  are  potentially 
useful  for  "rough"  engineering  eatimatt'S  of  sidelobe-to-eidelohe 
antenna  coupling. 

The  simplified  GTD  coupling  analysis  can  be  illustrated  by 
considering  the  coplacisr  arrangement  of  two  circular-aperture  horn-fed 
paraboloidal  reflector  antennas  shown  in  Figure  3- 17.  The  various 
symbols  shown  on  the  figure  are  defined  cs 

d^  aperture  diameter  of  Antenna  1, 
d2  *  aperture  diameter  of  Antenna  2, 
i j  *  length  of  OP, 

^2  “  length  of  OQj, 

13  •  length  of  OQ2, 

*  length  of  QjP, 

^•5  “  length  of  Q2P, 

”  separation  between  antenna  centers. 
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17.  Sketch  depicting  the  linear  distances  and  the 

angles  used  in  the  GTD  antenna  coupling  analysis 


The  angles  6^^  02,  &2t  ^  open  positive  as  shown. 

Two  different  types  of  coupling  arc  depicted  in  Figure  3-17. 
These  are  (1)  direct  coupling  of  the  feed  horns  along  OP,  (2)  and 
single  diffraction  coupling  along  OQj-QiP  and  also  along  0Q2~Q2^- 
Thus,  there  are  three  contributions  to  the  total  coupling  between  the 
antennas . 

The  author  of  Reference  24  calculates  the  coupling  coefficient  for 
each  type  of  coupling  when  only  that  particular  coupling  mechanism 
exists.  The  coupling  coefficient  is  defined  as  the  ratic  of  the 
received  power  to  the  transmitted  power.  The  total  coupling 

coefficient  q  when  both  types  of  coupling  exist  (implying  three 
contributions  as  noted  in  the  preceeding  paragraph)  is  computed  as 


n  •  ni'‘-n2Tn3  +  2Vn^  cosC^^^)  (3-22) 

+  2Vn,n7  cos(c,-)  +  2Vn,n  cos(C  ) 

23  23  X  3  13 

where 

■  power  coupling  coefficient  for  direct  OP  coupling, 
n2  •  power  coupling  coefficient  for  single  diffraction 
OQi-QiP  coupling, 

>13  ■'  power  coupling  coefficient  for  single  diffraction 
0Q2-Q2^  coupling. 


and  where  %2>  ^23*  ^13  relative  phase  angles  between  the 

electric  fields  arriving  at  P  via  paths  indicated  by  the  subscripts. 

The  power  coupling  coefficients  rij,  fi2»  '^3  presented  in 

Reference  24  ssy  be  written  in  our  notation  as 

<“2^  ]  (3-23) 


, 


H^(ai)  • 
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(3-24) 


Ho  - 


1611 


j^Sj(Y^,0^,0^)  •  £±11^(62).  an*! 


\  f' 

f<  v»  ,  o 


12 


)  '  4 

1611 


r_  -*1^2 

|s^(e^,0.  .0^)*  H2(Y2)J  sin  (.  ^)  . 


(3-25) 


where 

A  ^  the  wavelength  in  free  space, 

Gj  “  the  peak  gain  of  feed  horn  1,  and 
G2  the  peak  gain  of  feed  horn  2. 

Hi  and  Ho  in  Equation  (3-23)  are  the  vectorial  horn  patterns  with 
aziro.^thal  and  elevation  vectorial  components  iji  and  0  .  The  vector  Sj 
and  the  vectorial  pattern  H2  appearing  in  Equation  (3-24)  are  defined 
with  respect  to  the  edge  diffraction  point  coordinate  system  erected  at 
Ql  with  tangent  to  the  edge  and  0^  perpendicular  to  the  edge. 

Similarly,  S2  and  H2  appearing  in  Equation  3-25  are  defined  with 
respect  to  the  edge  diffraction  point  coordinafe  system  erected  at  Q2 
with  ^2  tangent  to  the  edge  and  0^  perpendicular  to  the  edge. 

Sj  and  S2  are  then  written  as 


)  =  \  (0i,03}  +  Hg  (6j)  F^j(0i,0s)  0i 

^  “  (3-26) 


and 

S2(8i.».)  -  *2  *  “2 

where  the  notation  H^^ ^  indicates  the  component  of  Hi  with  respect  to 
the  coordinate  system  (^i,®i),  etc.  The  angles  of  incidence  Oi  and  the 
angle  of  scatter  %  appearing  in  the  edge  diffraction  function  F  are 
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an^ociated  with  an  «dge  point  are,  of  course,  specified  with  respect  to 
;;he  particular  edge  point  coordinates  shown  in  Figure  3-18. 
ant  are  conputed  as 


P<M(®i.®s)  - 


■V 

p.l 

cos 

e.  ^ 

sin 

H 

cos(e  )  -  cos(6,) 
s  i 


(3-26) 


Fei(0i,08)  - 


cos (6  )  -  cos(?,)  ' 

s  1 


(3-27) 


where  it  is  assumed  that  the  scattering  angle  of  interest  is  away  from 
the  shadow  boundary  defined  by  •  *®i»  This  restriction  can  be 
eliminated  via  the  use  of  iscdified  edge  diffraction  fomilas  as 
presented  in  References  [25,26].  However,  (3-26)  and  (3-27)  are 
satisfactory  for  our  purposes. 

fhis  simplified  version  of  the  very  powerful  GTD  analysis 
technique  can  be  sed  to  make  rapid  computer  estimates  of  out-of-band  as 
well  as  in-ba  .d  antenna  coupling  by  specifying  the  feed  horn  patterns 
cf  the  two  antennas  over  the  frequency  intervals  of  interest.  It  is 
rgain  noted  that  the  equations  are  applicable  only  to  sidelobo-to- 
sidelobe  coupling  situations.  Stochastic  equations  based  on  the 
foregoing  CTD  equstions  could  be  derived  expressed  in  terms  or  the  mode 
excitation  coefficients.  However,  the  FVS  and/or  the  StfS  techniques  are 
envisioned  as  the  more  promising  methods  for  stochastic  coupling 
analysis  for  arbitrary  situations  and,  consequently,  the  GTD  analysis 
was  not  pursued  f^^rther. 
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SECTION  IV 

SYSTEM  DEVICE  EFFECTS 


A,  Introduction 

The  t«o  aain  objeecivea  for  this  task  tiere  (I)  to  develop  a  method 
for  analysing  the  antenna  pattern  effects  of  higher  order  modes 
generated  and/or  propagated  by  common  feed  system  devices  and  (2)  to 
assess  the  feasibility  of  deducing  the  mode  excitations  for  a  radiating 
feed  system  from  analysis  of  measured  out-of-band  pattern  data.  The 
first  objective  has  been  achieved  via  the  theory  and  equations 
developed  during  the  course  of  the  research  work  on  Task  I  and  Task  II 
preaented  in  Sections  It  and  III,  respectively.  The  results  of  the 
research  tiork  pertaining  to  the  second  objective  are  presented  and 
discussed  in  subsection  B  of  this  section. 

The  timely  completion  of  the  required  research  work  for  this  task 
permitted  a  brief  digression  into  another  important  aspect  of  the  out- 
of-band  response  of  system  devices.  In  particuler,  equations  were 
sought  to  describe  the  influence  of  particular  waveguie'e  devices  on  the 
excitation  and  subsequent  propagation  of  out-of-band  waveguide  modes. 
Three  waveguide  devices  were  considered  in  this  exploratory  study, 
nasKly  (1)  coax-to-waveguide  adapters,  (2)  radial  bends,  and  (3) 
ferrite  phase  shifters.  Program  constraints  precluded  a  detailed 
quantitative  analysis  of  all  three  kinds  of  devices.  Accordingly,  the 
research  efforts  were  concentrated  mainly  on  Che  theoretical  and 
numerical  analysis  of  a  coax-to-waveguide  adepter.  The  results  of  the 
exploratory  study  ere  suamurised  in  subsection  C. 

B.  Computatioo  of  Mode  Coefficients  Prom  Measured  Out-of-Bend 

Pattern  Pate 

The  mode  coefficients  for  the  propagating  modes  of  a  rsdiating 
feed  system  can  be  determined  from  e  knowledge  of  the  aperture 
tangential  electric  fields.  The  aperture  fields  can  be  determined  to 
good  approximation  from  Che  complex  far-field  electric  field  patterns 
of  the  feed  system. 

The  hypothetical  experimental  arrangement  consisting  of  an 
arbitrary  feed  system  terminating  in  a  large  metallic  "ground  plane"  is 
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phovn  in  Figure  4-1.  The  feed  system  is  assumed  to  consist  of  (1)  an 
input  sect ion  which,  for  example,  might  be  a  coax-to-waveguide  adapter 
or  a  powar  coupling  iris,  (2)  passive  and/or  active  devices,  such  as  a 
radial  b'>nd  and  a  phase  shifter,  and  (3)  the  radiative  element,  which 
is  an  unflared  straight  waveguide  section.  Note  that  the  guide  axis  is 
along  Che  x  coordinate  and  that  tne  y  and  z  coordinates  are  the 
transverse  coordinates.  This  orientation  of  cooidinate  axes  is  a 
convenient  one  for  the  probleo  at  hand. 

A  measuresKot  is  performed  to  determine  the  far-field  complex 
V  .torial  electric  field  in  the  forward  hemisphere.  The  aiessuremenC  is 
probably  most  efficiently  accomplished  via  Che  planar  near-field 
measurcoenc  technique.  The  far-field  electric  field  is  then  accurately 
computed  via  established  PFT  processing  of  the  near-field  data  [25]. 

The  transverse  electric  fields  B(y,z)  in  the  aperture  of  the 
radiating  system  may  be  computed  accurately  as 

E(y,B)  *  Tij  i^q  f-pq  ®*p[~jthypq  y  +  k^qt)  1  (4—1) 

where  che  coefficient  Epq  is  the  coaqtlex  far-tield  electric  field  in 
Che  direction  defined  by  (ky  *  ^ypq>  ^zq^  where 


■ypq  “  s  in(  Qq)sin\  4^)  , 

(4-2) 

*q  -  koCOsOq),  and 

(4-3) 

o  •  2TT/X 

(4-4) 

where  X  is  Che  free  space  wavelength. 

The  mode  coefficients  are  determined  by  solving  the  simultaneous 
equations  obtained  from  Che  continuity  of  Che  tangential  compoaents  of 
the  E  and  S  fields  at  the  aperture.  The  appropriate  equations  at  Che 
aperture  (x*0)  may  be  written  as 
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Sketch  depicting  a  rec-.angular  waveguide  transmission  line 
terminating  In  a  large  "ground  plane"  and  showing  selected 
geometrical  variables. 


(4-6) 
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where , 


a  >  the  coBplex  excitation  coefficient  for  the 
allowed  waveguide  «ode, 

b^  *  the  coaplex  reflection  coefficient  for  the 
^th  waveguide  node, 

_(+) 

y  (y.z)*  the  traoawerse-to-x  vectorial  waveguide 
electric  field  function  for  propagation 
in  the  (*)  x-direction. 
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the  effective  wave  adaaittance  dyadic  for 
the  radiated  fields,  Fpq  ~  1^,  sin(eq) 
cos  6p/utJ  ,  and 
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the  wave  adaittance  dyadic  for  the  Mveguide 
field  functions. 
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where  Fpjj  «  [kQ8in(8q}cos(6p)  ]/«p  - 
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The  waveguide  electric  field  functions  are  assuaed  to  be  arranged 
in  ascending  order  of  cut-off  frequency.  The  electric  fields  for  the 
TE  and  TM  aodes  that  have  the  atme  cutoff  frequency  require  special 
handling  in  order  to  preserve  a  desirable  orthogonality  prorerty  of  the 
node  fields.  In  particular,  if  (y>x)  is  :i  TE  aode  field  with  cut¬ 

off  wavelength  and  (yts)  is  the  TM  aK>de  field  with  the  aaate  cut¬ 
off  wavelength,  the  waveguide  field  functions  (y,*)  and  are 

teken  to  be 

4  _+ 

i~(y.z)  -  U^(y,z),  and  (4-10) 

+  +  + 

‘  *  V‘(y.z),  (4-U:. 

where 


a 


K 


j(/tu; (>.«))*• 


+ 

V^(y, z)dyd2 
~ 


U^(y,2)dydz 


(4-12) 


is  set  equal  to  aero  for  any  pair  of  nodes  that  have  different  cut¬ 
off  wavelengths.  The  integrals  are  evaluated  over  the  guide  cross 
section  • 

The  a-polariaed  and  y-polarized  electric  fields  for  the  TE 
wsveguide  iBodes  traveling  in  the  positive  x-direction  are 
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luICy.*)] 


*  [{~f^  +  (— )^1 

D  9 


sin  (y  +  f)  1 


co8[^  (a  y)l  exp  (^-13) 
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D  a 
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(4-14) 


where 


Vk^ 


(5,,  2 

a 


b"' 


aad 


(4-35) 


where  the  n  and  m  indices  are  chosen  based  on  the  node  ordering  schetse 
previously  dicussed.  Siailarly,  the  s-'polarized  and  y-polarized 
electric  fields  for  the  TM  nodes  are 


tv“(y,z)ljj 


!(“)*  -  (T)'l 

D  a 


sin  [“  (y  +  |)  ] 


cos  (2  +  j)]  exp  and 


(4-16) 


[v^(y.z)]^ 


.  -in'*,  2 

I'T'  ♦ 


A 


cos  [ 


ns 


sin  1“  (x  +  1^)]  exp  [-j8^*l 


(4-17) 


It  should  be  noted  here  thet  node  seta  other  than  the  node  sets  ?B  and 
TM  to  the  longitudinal  axis  can  be  employed.  However,  this  is 
frequently  a  convenient  set  of  nodes  and  serves  to  illustrate  the  basic 
theory  and  equations. 
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the  solution  for  tue  node  coeff icient.s_  is  readily  achieved  by  (1) 
scalarly  vultiplying  Equation  {4~6)  by  integratiug  both 

sides  of  the  equation  over  the  guide  cross  section :  (21  scalarty 

atultiplying  Equation  (4-5)  by  (f  *  ***^  integrating  both  sides 

of  the  equation  over  the  guide  cross  section,  making  use  of  the 
orthogonality  of  the  nodi  fields  to  retain  only  the  <  *  id  terms,  and 
solving  the  simultaneous  equations  in  the  unknown  a  and  b  for  a 

K  K  K 

The  solution  for  a  ao  obtained  ia 


a 

K 


I  z  i 

p  q 


a  b 


2  2 


ir  +  r  ) 
_ K _ ES— 


E 

_ E3..- 

•  f^(y,*)drdy 


(4-18) 


where 


a  Jb 

-  J  ,  v,-+  *  (4-19) 

“b. 

y  2 

is  the  complex  conjugate  of  the  vectorial  far-field  electric  field 
pattern  radiated  by  the  waveguide  mode.  All  of  the  quantities  on 
Che  ri^t-hand  side  of  Equationa  (4-18)  and  (4-19)  ace  known  from 
experiment  or  theory  and,  hence,  the  coefficients  may  be  computed. 

C.  Sumnary  of  Out-of->and  Wavcgiide  Device  Characteristica 

The  emphasis  in  this  brief  study  was  on  oti'aining  a  qualitative 
understanding  of  the  basic  out-of-band  mode  ex;citations  and  propagation 
characteristics  of  (1)  a  coax-to-vsveguide  adapter,  (2)  radial  bends, 
and  (3)  ferrite  phase  shifter.  A  detailed  quantitative  study  of  each 
device  wa»  outside  the  scope  of  these  exploratory  investigations  due  Co 
the  complexity  of  the  solutions  obtained  for  mulcimoding  waveguides. 
Accordingly,  attention  was  focused  primarily  on  the  coax-to-waveguide 
adapter.  The  key  analytical  results  Cor  the  radial  bends  and  the 
ferrite  phase  shifter  are  first  briefly  aummarixed  in  the  following 
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paragraphs.  rhe  salient  ;inalyeical  nuBerical  results  for  the  coax-to» 
waveguide  adapter  ace  then  presented  and  discussed. 

1 .  Waveguide  Radial  Bends 

The  propagation  of  electraujgnetic  energy  through  radial  **E- 
plane”  and  "H-plane"  waveguide  bends  hu:  )k  er  studied  by  a  large  number 
of  authors  during  the  last  45  years.  A  partial  list  of  references  is 
included  in  Section  VII  as  References  26  throirgh  34.  Out-of-band  as 
veil  as  in-band  propagation  has  been  addressed.  Rigorous  solutions  for 
the  electric  fields  in  beads  are  presented  in  References  26,  31,  32, 

and  33  in  tersis  of  infinite  series  of  Ressel  functions.  The  matrix 

solution  due  to  Rice  [30]  is  also  capable  of  providing  rigorous 

results.  Rigorous  solutions  for  the  generation  of  higher  order  modes 

ct  the  iunccion  of  straight  and  curved  portions  of  a*  waveguide  run  have 
else  been  aerived  [32].  The  rigorous  solutions  are  generally  valid  for 
"sharp"  bends  as  well  as  for  "gradual"  bends.  A  sharp  bend  it  defined 
herein  as  one  for  which  AR/R  i  1  and  a  gradual  bend  is  onii  for  which 
AR/R  <<  1,  where  R  is  the  bend  radius  and  A  K  is  the  thickness  of  the 
bend  alo\vr.  ? ,  as  per  Figure  4-2. 

Thr  approximate  analysis  of  a  gradual  H-plrne  band  by  Jouget 
[27, 2C]  furnishes  valuable  physical  insight  as  well  as  useful 
approximate  equations  for  the  intermodsl  coupling  coefficients  at  the 
junction  of  the  straight  and  curved  guides.  The  corresponding  analysis 
for  the  round  waveguides  are  also  presented  by  Jouget  to  a  seps: ate 
paper  [29]-  Expressions  for  the  mode  cutoff  frequencies  and  the  mode 
phase  velocities  in  the  curved  portions  of  the  wav'^guide  are  also 
derived.  The  corresponding  analysis  for  an  E-plane  bend  is  presented 
to  second  order  by  Lewin  [34].  Good  engineering  approximations  for  the 
mode  conversion  coefficients  at  the  entrance  and/or  exit  junctions  are 
possible  even  ‘:or  aod:!rately  sharp  tenda  where  AR/R  t  c.X. 

2.  Ferrite  Phase  Shifter 

The  :wiin  difficulty  in  the  analysis  of  a  waveguide  containing 
a  ferrite  pha'^e  shifter  is  in  obtaining  a  valid  repreaentation  for  Che 
electric  field  or  Che  magnetic  field  in  Che  portion  of  the  guide 
containing  Che  ferrite  slab.  Once  these  cbaracCerisCic  fields  are 
obtained,  the  usual  proceduie  of  requiring  continuity  of  the  fields  at 
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geometrical  varlablee. 


the  entrance  an4  exit  *aces  of  the  slab  can  be  invoked.  This  leads,  ot 
course,  to  a  set  of  sisultaneous  equations  that  are  readily  solved  to 
provide  reflection  and  transmission  coefficients  as  functions  of  the 
phase  shifter  control  setting  for  a  specified  incident  eode.  For 
exaaple,  the  fields  for  a  rectangular  guide  containing  a  ferrite  slab 
extending  fully  across  one  diaension  of  the  guide,  as  per  Figure  4~3, 
can  be  solved  exsctly.  Houever,  the  exact  solution  for  slab  geoaetry 
shown  in  Figure  4-4  is  too  coaplex  to  be  useful  and,  hence,  approximate 
techniques  are  sought. 

The  ferrite  phase  shifter  is  generally  a  lossy,  anisotropic  aediua 
and,  consequently,  the  application  of  perturbational  or  variational 
techniques  is  aore  coaplex  chan  for  typical  waveguide  probleas  [  11,33]. 
Nonetheless,  the  application  of  variational  techn'ques  to  the  ferrite 
phase  shifter  problea  is  potentially  very  useful.  The  equations 
presented  by  Hoshen  [36]  appear  to  be  particularly  well-suited  for  this 
problea.  Moshen  derives  a  stationary  quantity  P  that  ia  valid  for  a 
lossy,  anisotropic,  inhoaogeneous  aediiai  which,  for  a  source  free 
region,  is 


F  -  [V  X  E.H]  +  [V  X  H,E] 
s  s 


+  ju  ([B.H]^  -  (O.EJ  ) 

S  5 


(4-20) 


+  lE.H]^ 

where  quantities  of  Che  fora  [A,!]^  and  [  4,8]^  denote  Che  "inner 
products" 


ds  and 


(4-21) 


ilc-  / 


A  •  (n  X  B  do. 


(4-22) 


Figure  4-3.  Sketch  depleting  «  section  of  rectangular  waveguide  containing 
lerrlte  waveguide  slab  whose  height  equals  the  waveguide  heigh: 


Sketch  depicting  e  section  of  rectangular  waveguide  containing 
a  territe  slab  whose  two  transverse  dimensions  are  less  than 
the  corresponding  dlmenalons  of  the  waveguide. 


rei^pectively.  The  synbola  E,  H,  and  B  denote  the  electric  field 
vector,  the  magnetic  field  intensity  vector,  and  the  magnetic  field 
vector,  consecutively.  The  inner  product  [  A,B]g  is  evaluated  over  the 
guide  cross  section.  The  inner  product  is  evaluated  around  the 

periphery  of  Che  guide  cross  section.  The  suoerscript  syisbol  t 
denotes  the  adjoint  of  th<*  field  quantity  to  which  it  is  appended. 

A  variational  expression  for  the  propagation  constant  in  the 
ferrite  pha^e  shifter  can  be  obtained  from  Equation  (4-20).  Tue  trial 
field  E(x,y,s)  may  be  conveniently  chosen  as  the  superposition  of  mode 
fields  f^(x,y,*)  for  the  empty  waveguide,  to  wit: 


E(x,y,*)  - 


(4-23) 


The  irial  U  and  B  vectors  are  obtained  from  Equation  (4-23)  via  the 
Maxwell  cuil  relation  and  the  constitutive  relation  B  ’*  u  *  B  ,  idiere 
P  is  the  permeability  tensor  [34]  for  the  ferrite  slab.  The  expansion 
coefficientc  which  extremixe  the  propagation  coefficient  are  found 
by  numerically  solving  the  set  of  simultaneous  equations  obtained  from 
the  variational  equations.  Thus,  the  electric  field  E  given  by 
Equation  (4-23)  is  determined  and  can  thence  be  employed  in  the 
analysis  and  computation  of  reflection  and  transmiiision  coefficients 
for  the  ferrite  section. 


3.  Study  ot  the  Coax-to-Waveguide  Adapter  Device 

The  goal  of  this  study  is  to  obtain  a  qualitative 
understanding  of  the  higher  order  mode  excitation  and  propagation 
characteristics  of  the  coax-to-waveguide  adapter.  The  results  of  this 
study  provide  useful  insight  into  the  basic  out-of-band  characteristics 
of  a  coax-to-waveguidc  adapter. 


The  solution  for  the  mode  excitstion  coefficients  for  the  coax-to- 
waveg’iide  adapter  shown  in  Figure  4-5  can  be  obtained  with  the  aid  of 
equations  given  in  Collin  [34].  Note  that  y  and  x  are  transverse 
coordinates  sad  the  longitudinal  axis  of  the  guide  is  along  z.  In 
particular,  the  dyadic  Green  function  C(x' ,y* ,z' |x,y (Z)  derived  by 
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Sketch  depicting  a  side  view  of  a  coax-to -waveguide  adapter  and 
showing  selected  geometrical  variables. 


Collin  for  an  inf initesiaal  electric  dipole,  located  at  (x*,y',*')  in 
3  rectangular  viaveguide  is  applicable  to  both  in-band  and  oot-of-band 
waveguides.  The  aagnetic  vector  potential  A  •  Ayy  produced  by  an 
assuated  current  density  J  “  y  Josin  [k,>(d-y')]  on  the  probe  and  is 
derived  by  integrating  the  product  6  *  J  along  the  probe  surface  Sq,  to 
»it: 

A<x,y,z)  “  S  ^^*' *^' ’**  ^*’^**^  "  (4-24) 

Expressions  for  E(x,y,z)  and  B(x,y,z)  are  then  readily  obtained  via  the 
well  known  operations 

E(x,y,z)  *  -jw  A(x,y,z)  ^(7*4)  (4-25) 

jwc  y 
o  o 

and 


H(x,y,r)  -  —  (VxA)  (4-26) 

The  electric  field  l(x,y,z)  obtained  frosi  Equation  (4-25)  can  be 
%rrittea  as  an  expansion. 


E(x,y,z} 


I 

n 


Z 

m 


b 

am 


Xn„(x.y.t) 


(4-27 ) 


where  ■  -  0, i ,2,3,4, ... .  and  n  ■  odd  positive  integer  n  *  I, 3, 5, 7,  ...., 
the  selection  of  the  values  of  index  n  being  imposed  by  the  requireaent 
that  the  vector  potential  X(x,y,z)^  0.  The  expansion  functions  Xm  are 
the  TMj^  waveguide  mode  functions  which  have  no  coaponent  of  H  along 
the  transverse  y  direction.  (Thus,  one  can  anticipate  that 

longitudinal  power  flow  is  associated  only  with  the  y-polarized 

electric  fields  of  the  TMj^  aedes  since  *  0  ) .  The  coefficients 

are 
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exp[-j6  z]  -  exp[-j8  (z-M)] 

OB _ HA 


JB, 


na 


[ 


It  [cos (k  d)  -  cofci(i»«  — )] 

2r.  sinl^lj  - 2.__ - - - b - ^  ^ 


2  I  a 


(E)2  - 

b  o 


M.  • 

4 


where  *■  i  for  *  *  0 
and  ^o»  *  2  otherwise. 


The  vector  function  XQB(x,yt£)  is  expressed  in  tersu  of  its 
as 


•  sin  ■*•  y)]}x 


2y 


1 


o  0 


•  sin  (x  +  f^)]  cos  (y  +  ^)])j. 


Xnae*  “ 


(*  +  f>J 

D  niB  a  z 


cos  (y  +  y)] }*. 


(4-28) 


components 


(4-29) 


(4-30  ) 


(4-31  ) 
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The  electric  field  E(x,y,s)  can  be  expressed  with  the  aid  of  Equation 
(4-25)  aa 


,  ,2.  .  . 

'Ii3f  •  +  ^  *  5J57  '  > 


(4-32) 


H(x,y,z}  via  Equation  (4-26)  as 


B(x,y,z)  -  —  *  -  3^  ^ :) 

o 


(4-33) 


The  average  power  density  transported  along  the  z  direction  is 


A  ^  i  vA  ■ 

S*  -  S  •  2  -  ^  Re  E  .  (-^)  1 


y  3i 


(4-34) 


where  S  la  the  pointing  vector  defined  as 


S  »  y  Re (1x8*) 


(4-35) 


d*  the  cosiplex  conjugate  Magnetic  field  vector  and  s  the  unit  vector  in 
the  direction  of  propagation.  Utilising  Equations  (4-24  to  4-36),  a 
closed  form  expression  for  the  power  propagated  in  a  coax-to- 
waveguide  adapter  can  be  derived  by  integrating  Equation  (4-35)  over 
the  guide  cross  section  Co  obtain 
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where  Re  aeans  "real  part  of"  and  triiere  all  other  ayabola  have  been 
defined  previously,  the  total  power  propagating  along  the  waveguide  is 
partitioned  asong  the  different  possible  wodes  that  have  n  >  0  for 

a  particular  operating  frequency.  Calculations  of  the  power 
distribution  aaong  the  different  higher  order  Modes  were  made  for 
standard  S-band  coax-to-waveguide  adapters  for  frequencies  ranging  from 
3.0  to  10.0  CHz.  The  values  of  the  various  physical  parameters  of  the 
adapter  and  waveguide  for  the  numerical  study  are  as  follows: 

a  »  7.112  cm 
b  >  3.302  cm 
d  *  1.9  cm 
z  »  3.5  cm 
I  »  2.4  cm 

The  ca'iculations  were  performed  on  the  asstaption  of  the  total  power 
radiated  to  be  equal  to  1.0  vatt  corresponding  to  a  dipole  resistance  R 
=  500  ohms  and  a  current  «  43  milliamps. 

It  should  be  noted  again  that  the  power  flow  down  the  guide  is 
associated  entirely  with  the  y-polarised  component  of  the  electric 
field  of  the  excited  modes.  The  x-component  of  the  electric  field  does 
not  contribute  to  the  longitudinal  power  flow  since  the  y  componet  of 
the  magnetic  field  is  aero  for  the  TMA  modes. 
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Results  of  the  computer  calculations  are  tabulated  for  selected 
frequencies  in  Table  4-1.  In  particular,  the  modal  power  flow  for  each 
excited  TM^  modes  is  tabulated  as  the  percentage  of  the  total  power 
flowing  in  all  of  the  excited  modes^for  the  in-band  frequency  of  3.U  GHz 
and  for  the  out-of-band  frequencies  of  6.0  GHz,  7.5  GHz,  9.0  GHz,  9.5 
GHz,  and  10.0  GHz.  Inspection  of  the  tabulated  data  shows  how  the  total 
power  flow  is  partitioned  among  the  excited  modes  for  each  frequency. 
The  dispersive  character  of  the  waveguide  over  the  out-of-band  frequency 
range  from  6.0  GHz  to  10.0  GHz  is  clearly  evident  in  the  data.  The 
frequency  sensitivity  of  the  modal  power  flow  distributions  for  the 
frequencies  of  6.0  GHz,  7.5  GHz,  9.0  GHz,  and  9.5  GHz  is  partially  due 
to  the  fact  that  the  total  power  flow  is  partitioned  among  2  modes,  3 
modes,  4  modes  and  5  modes,  consecutively.  However,  considerable 
frequency  sensitivity  is  also  evidenced  for  the  out-of-band  frequencies 
of  9.5  GHz  and  10.0  GHz,  both  of  which  have  the  total  power  flow 
distributed  among  the  same  5  TM^  modes. 

It  should  be  noted  that  the  analysis  of  a  coax-to-waveguide 
adapter  having  a  non- filamentary  current  probe  can  be  obtained  via  the 
equations  presented  in  Reference  35  and  herein.  The  Green  function  for 
empty  rectangular  waveguide  presented  in  Chapter  7  of  Reference  35  must 
be  numerically  integrated  over  the  surface  of  the  current  probe  that  is 
being  modeled.  It  would  be  particularly  beneficial  to  analyze  the  out- 
of-band  behavior  of  a  coax-to-waveguide  adapter  having  the  shape  of  an 
elongated  "teardrop"  found  in  many  commercial  adapters.  Adapters  which 
have  a  concentric  dielectric  sleeve  enclosing  the  current  probe  require 
the  use  of  a  more  complicated  Green  function.  It  would  also  be  of 
considerable  practical  value  to  numerically  analyze  the  out-of-band 
characteristics  of  the  dielectrically-sleeved  adapters. 
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TABLE  4-1 


MODAL  POVER  FLOW  FOR  TBE  MODES  FOR  THE 

INDICATED  FREQUENCIES  FOR  WR-284  VAV^IOE 
EXCIiED  Bf  AR  IDEALIZED  COAX-TO-HAVBGOIDE  ADAPTER 
HAVING  A  FILAMEHTARY  CURRENT  FROBE« 


Frequency 

Modes 

Modal  Power  Flow 

(GHz) 

(n,iB  indices) 

(X  of  Total) 

3.0 

10 

100.0 

6.0 

10 

6.96 

11 

93.04 

7.5 

10 

11.99 

11 

7.00 

30 

81.01 

9.0 

10 

47.05 

11 

13.62 

30 

0.45 

31 

38.88 

9.5 

10 

47.83 

11 

17.86 

30 

10,64 

31 

6.40 

12 

17.27 

10.0 

10 

32.86 

11 

11.51 

30 

20.33 

31 

0.00 

12 

35,30 

waveguide  aodes  are  defined  with  respect  to  the  coordinate 
syTitem  depicted  in  Figure  4-5. 
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SECTION  V 
SITE  EFFECTS 


A.  Introduction 

It  is  well  known  that  the  performance  of  an  antenna  can  be 

significantly  degraded  due  to  scattering  of  electromagnetic  energy  from 
metallic  or  dielectric  structures  located  within  the  radiating  near¬ 

field  of  the  antenna.  Several  techniques  are  available  for  analyzing 
the  effects  of  near-field  obstacles  for  a  specified  operating 
frequency.  Prominent  among  these  techniques  are  aperture  integration, 
the  Geometrical  Theory  of  Diffraction  (GTD),  and  the  Plane  Wave 

Spectrum  Scattering  (PWSS)  analysis.  Although  any  of  these  techniques 
could  be  used  for  out-of-band  frequency  analysis,  they  have  generally 
been  applied  to  the  analysis  of  in-band  effects  at  a  few  selected 

frequencies  of  interest. 

The  analysis  of  the  distortion  of  the  envelope  of  a  narrow  time 
pulse  requires  a  knowledge  of  the  frequency  domain  antenna/obstacle 
response  over  a  frequency  inteval  exceeding  the  bandwidth  of  the  pulse. 
Alternatively,  a  direct  time  domain  analysis  can  be  formulated  and 
applied  to  near-field  scattering  problems,  as  per  Reference  [37  ]  . 
However,  a  frequency  domain  analysis  appears  to  be  more  flexible  and 
appropriate  for  the  wide  variety  of  antenna/scattering  situations  that 
are  encountered  at  Army  directive  antenna  installations.  Moreover, 
experience  has  shown  that  the  PWSS  method  of  antenna/obstacle  analysis 
is  particularly  well  suited  for  analyzing  near-field  obstacle  effects. 
Accordingly,  the  frequency  domain  approach  employing  the  PWSS 
scattering  technique  has  been  used  to  study  the  effects  of  near-field 
obstacles  on  the  wideband  CW  and  pulsed  response  of  directive  antennas 
for  both  in-band  and  out-of-band  frequencies.  The  basic  theory  and 
equations  are  presented  and  discussed  in  subsection  B,  and  the  results 
of  a  numerical  study  of  wideband  site  effects  are  summarized  in 
subsection  C. 

B.  PWSS  Formulation 

Tne  monochromatic  PWSS  analysis  is  readily  adapted  to  analyze 
near-field  obstacle  effects  for  wideband  CW  or  pulsed  systems.  The 
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theory  and  equations  for  the  deteministie  MOaochr<matic  PVSS  are 
outlined  first. 

1 .  Plane  Wave  Spectnai  Analyaie  Coccept 

!nie  basic  theory  of  the  Plane  Wave  Spectnas  Scattering  Matrix 
analysis  is  presented  in  Keferenc'^s  22  and  38  .  A  brief ^  sonewhat 
heuristic  synopsis  it  presented  here.  Figure  5-1  illustrates  the  basic 
concept  and  the  antenna/obstacle  geoaetry.  Conceptually,  the  antenna 
Plane  Wave  Spectrus  (JVS)  referred  to  the  antenna  aperture  propagates 
in  the  near- fie  Id  and  is  incident  upon  an  obstacle.  Each  cossponent 
plane  wave  in  the  antenna  spectrun  excites  a  conplete  spectrun  of 
obstacle  scattered  plane  waves.  The  total  scattered  apcctrua  is  thus 
the  superposition  of  the  scattered  plane  wave  spectra  due  to  the 
incident  antenna  spectrun. 


If  the  antenna  aperture  is  located  in  the  x  “  0  plane,  the  antenna 
FVS  A(ky,k2}  can  be  expressed  as, 


j(k  ys-k  z) 

e  ^  dydz,  (5-1) 


where  k^  •  Vg  sinG  cos  (J)  ,  ky  »  ko  sin9  sin  <|i  ,  k^  ■  k^  cos  8  ,  ■ 

2’t/^,  E[(0,y,r)  *  yEy(0,y,z)  +  i»  tbe  electric  field 

transverse  to  x  and  an  e^i"^  tine  dependence  is  suppressed. 

Equation  (5-1)  expresses  the  transverse  antenna  specCrua  A(ky,k2) 
as  th»>  Fourier  trans format ion  of  the  transverse  aperture  electric  field 
•  ''ia  the  inverse  transforaation  the  tranaverss  to  x  electric  field 
in  the  charge-free  re<,,ion  x  >0  can  be  expressed  as  111, 35  } 


dk  dk  . 

y  * 


(5-2) 


where  k  »  xk,{  v  yky  +  zk^  and  r  “  fr  ■  xx  +  yy  +  zz. 

In  the  far-ficld  where  x  0,  an  asymptotic  expression  can  be 
derived  by  the  method  oi  steepest  descent  as, 

E^(r.9.d)  •  r  * 


i:o 


where  the  wsvenuabere  ky^  and  are  the  values  of  ky  and  k^  at  the 

saddle  point  of  the  integral  of  E<;uation  (5-2). 

The  far-zone  scattered  field  due  to  a  tingle  conponent  of  the 
plane  wave  spectrua  incident  on  Che  obstacle  shown  in  Figure  5—1  can  be 
expressed  as 

m 

E®(r\e.*)  -A  {ty.k^;kj.k^)  •  A^(kJ.k^)  .  (5-4) 

where  A^(ky,V*)  -  A(ky,kj)  e"iki  -8,  f  denotes  the  obstacle  location 
and  is  the  plane  wave  to  spherical  wave  scattering  dyad  for  the 

obstacle  with  phase  rifercnce  at  the  obstacle  location.  The  incident 
wavenuaiber  k^  is  given  by 


-  i  /k^2  .  _  (kj)*  +  +  Jkj  .  (5-5) 

and  the  term  ‘  ^  shifts  the  antenna  spectnsa  to  the  obstacle 

reference  center  located  at  R. 

For  our  purpose,  it  is  convenient  to  express  Equation  (5-4)  as 
i'Cky.k^)  -  SCky,k^;kJ,k^  *  i^(ky,kj)  .  (5-6) 

where  by  using  Che  asyaptoCic  relation  of  Equation  (5-3) 


-  fsk-^ 


(k  ,k  jk^.kS 

7  •  7  * 


(5-7) 


is  Che  plane  wave  spectna  scattering  dyad  that  specifies  the  obstacle 
scattered  plf-ne  wave  spectrua  for  the  case  of  excitation  by  an  incident 
plans  wave  traveling  in  the  direction  defined  by  (ky^k^). 

The  total  obstacle  scattered  plane  wave  spectrum  can  now  be 
expressed  as  a  superposition  of  the  scattered  spectra  due  to  each 
component  of  the  incident  antenna  spectrum  as, 
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A*(k  ,k  )  - 
y  * 


L  Si- 

f  f  •  <5-8) 

“‘‘.L  -Sl 


where  k^^  *t^  k^,  ere  the  vevcouaber  Imite  of  the  iocident  spectrua. 

In  practice,  it  is  aore  convenient  to  evaluate  Equation  (5~8)  in 
the  angle  (e  ,<» )  space  since  aost  computer  programs  for  scattering 
analysis  employ  (6,4)  coordinates.  This  transformation  results  in, 


■^•fL  TL 


•{6. 4)  -  -k^^  f  f  •  A(8^4^) 


^  2-i 

Bin  6  cos4  dd  d4  , 


(5-9) 


%dtere  6^  and  are  the  angular  limits  on  the  incident  spectrua  space. 

It  is  also  convenient  to  express  the  scattering  dyad  in  rersis  of 
spherical  coord inants  as 


•  1  i 

/  ♦♦ 

♦e 

*48 

S(e,4;e^,4^)  - 

i  -- 

»e 

*88 

(5-10) 


where  is  the  scattering  function  for  the  <4-po}arized  component  of 

Che  scattered  field  for  a  8 -polarized  incident  plane  wave  and  so 
forth.  It  is  apparent  that  this  form  fits  nicely  with  scattering 
analyses  expressed  in  terms  of  tl'e  conventional  TE  and  TM  incident 
plane  wave  Bx>des,  and  the  scattert.d  spectra  for  Che  obstacle  are 
readily  expressed  in  terns  of  orthogonal  (8,4)  components  [  11  }. 

The  resulting  total  spectrum  is  given  by  Che  sum  of  the  antenna 
and  scattered  spectra  in  the  (8,4)  space  as. 
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-  «lAg<fl.*)+Aj(e.4)]  +  >tA^(e,4)+Aj(S.^)],  (5-11) 


and  Che  coopoeite  far~zone  pattern  can  be  readily  obtained  fros  the 
total  spectrum. 

The  specCrtOi  ACkyikf)  which  is  derived  via  the  Fourier 
transformation  of  the  transverse  to  x  electric  field  E(0,y,x)  is  tersted 
the  "transverse  vectorial  spectrum"  [21  ]  .  Since  fai-eone  radiation 
patterns  are  conventionally  expressed  in  terms  of  the  angular  (  9,  ij) 
coordinates,  and  it  is  sometiws  convenient  Co  evaluate  the  scattering 
matrix  for  (6,#)  coordinates,  it  is  often  necessary  to  convert  the 
transverse  vectorial  spectrum  expressed  in  the  (x,y,x)  cartesian 
coordinates  to  Ch£  angular  coordinates.  To  accostplish  this 

transformation,  we  define  after  Kerns [21  ],  Che  orthogonal  unit  vectors 


(5-12) 


(5-13) 


Kote  Chat  the  unit  vectors  Kj  and  K.2  are  respectively  "radial"  and 
''tangential"  unit  vectors  associated  with  polar  coordinates  in  the  (y,x) 
plane.  The  individual  plane  wave  components  of  the  spectrum  propagate 
in  the  plane  defined  by  Kj  and  the  x  axis.  Thus  an  electric  field 
aligned  with  Kj  is  parallel  to  the  plane  of  incidence  (Etm 
Transverse  Magnetic)  and  an  electric  field  aligned  with  K2  is 
perpendicular  Co  the  plane  of  incidence  vKie  or  Transverse  Electric). 

The  scalar  spectral  density  functions  corresponding  to  Che  basis 
vectors  K|  and  K2  outgoing  vectorial  spectna  ere  defined  as 


bd.K)  -  ACK)  , 

b(2,K>  -  *2  *  * 

i2i. 


C5-14) 


where 


K  “  yky  +  zk* 


Finally,  Che  complete  vectorial  Rpectrxxm  a(R}  i^ich  includes  all 
of  the  cartesian  (x,y,a)  components  of  the  outgoing  spectrum  is  given 
by  the  relation 

a(K)  -  b(l,K)[Ki  +  X  +  b(2,K)  K2  ,  (5-15) 

'x 

where 


The  completa  cartesian  baaed  vectorial  spectrum  a(R)  can  then  be 
converted  to  Che  angular  (r,^,$)  basis  vectors  via  the  relations 


X  ■  r  sine  cos^i  +  0  coa6  co3<j)  -  ^  sin^  , 

y  »  r  sine  sin^  +  0  cose  sin^  ♦  0  cos^i  ,  and  (5-16) 

z  *  r  CO80  -  0  sin6 

The  foregoing  analysis  which  was  described  for  a  single  mast 
obstacle  is  applicable  to  any  number  of  obstacles.  It  is  only 
necessary  to  specify  Che  scattering  dyad  S  for  the  obstacle  cluster 
and  the  vector  R  from  the  antenna  to  a  designated  phase  reference 
center. 

The  PWS  dyad  S  for  the  obstacle  or  group  of  obitacles  may  be 
specified  fiom  theory  or  experiment.  The  FWS  dyad  can  be  constructed 
^rom  GTC,  physical  optics  or  the  Method  of  Moments  solutions.  Tho  FWS 
dyad  for  amlciple  obstacles  must  account  fer  multiple  scattering  among 
tne  obstacles  in  order  to  achieve  accurate  results  for  closely  spaced 
scatcerers  [39  ]  .  However,  accurate  approximate  equations  are 
available  only  for  collecciona  of  parallel  cylinders  cr  spheres  of 
small  diameter.  Generalized  practical  techniques  for  constructing  the 
PWS  dyad  for  closely  spaced  obstacles  of  arbitrary  shapes  are  not 
available  at  the  present  time.  Accordingly,  this  remains  as  one  of  the 
most  important  problem  areas  in  the  analysis  of  site  effects  for  many 
real-world  antenna  installations. 

The  PWS  technique  is  readily  adapted  for  wideband  antenna  analysis. 
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The  extension  to  wideband  stochastic  in-band  and  out-of-band  situations 
is  presented  and  discussed  in  the  following  parasraphe. 

2.  Wideband  Analysis 

It  will  be  advantageous  to  %rrite  the  total  electric  field  for 
radian  frequency  ui  ■  2ir  f  as 


E(u).r.9,q)) 


ih(e.  <(>)exp[-i(.)(r/c)] 
cr 

+  exp(  ]  A®(m,0,ij>)] 


(5-171 


where 

E(  (ij.r.e  ,4>  )  ”  total  field, 

h  (  0,  <)>  )  ■  sin0  cosiji 

K  *  8  in  6  co8(i  x-t'sinBsini^  y  +  eosd  z 

c  *  the  speed  of  light  in  free  space 

where  the  superscripts  a  and  s  on  A  denote  the  antenna  and  scatterer, 
respectively.  The  scattered  antenna  spectrua  in  Equation  (5-17)  is 
referred  to  a  coordinate  origin  located  in  the  geoswtrical  center  of 
the  scatter.  The  exponential  factor  exp  [  jw  (K  '  E/C)]  "aoves"  the 
scattered  spectrun  to  location  of  the  origin  of  the  antenna  rotation 
coordinates.  The  quantity  (1/cr)  will  be  suppressed  in  the  reaaining 
equations,  and  tht.  total  field  will  be  written  as  E(  u,9  )  in 

recognition  that  the  factor  (i/cr)  is  being  suppressed. 

TTie  antenna  spectrum  function  A*(  u>  ,  6  ,  4*  )  for  an  ovensoded 
reflector  antenna  is  comprised  of  the  weighted  susi  of  spectrums  for  all 
of  the  propagating  modes  'r  the  antenna  feed  system.  Hence,  the 
resultant  antenna  spectrum  is 

A*(id  ,  0  , 4)  )  •  a^  (w  )7y  (  n),9  ,  4>  )  (5-18) 

where 

a  (ui  )  ■  the  complex  excitation  coefficient  for  the 
^  waveguide  mode  Tor  frequency  u,  and 


2 


(U)  ,0  -  the  antenna  spectrum  produced  by  the  waveguide 


mode. 


It  then  follows  that 


the  scattered  spectrum  function  A  (w  ,0  J 


A®(  0),  0,  <fr)  -  \ 


(5-19) 


where 


A®(ai.0.4>) 


.(iiyZ  y  y  {S 

\  -♦l 

•exp[-j  ^  [(sln(0)cos((>))]X  +  [sin(0’)sin(4.')  ]Y  +  lcos(e')]Z 


•sin^(e')cos(<|)')d(>’  d0’} 


(5-20) 


The  detemini.tic  pewer  d.n.ity  PC..,  e,^»)  .t  obeer.ation  angle. 
<  e,  ♦)  for  treqneney  »  i.  propoitlnn.l  to  E  (  « ,  6  .  ♦ )  •  EC  »,e  ,d  ) 
and  ia  written  explicitly  in  tern,  of  the  conplex  «nie  exettatton 

coefficients  a(w  )  as 


P(<*>,e,<|.)  =  a*.(a))  a^(u)) 


(5-21) 


where 


P  ,  (co,0,<t>) 


u)V(0.^)[A®?(w.e,«t')A®  («D,0.<l»)  +  A®,(w.0,^)A^(u».e.<l>) 


+  exp  [+jw  A®*(u>.0,«^)a“(u),0.<>) 

+  exp  l-ju.  A®(u..0.«)A®;(a..0.^)l  (5-22) 
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Th?  statistical  average  power  density  <  P(bi  ,6  )>i8  obtained  from 

Equi;tion  (5-21)  by  statistically  averaging  over  the  product  (m  ) 

of  the  complex  mode  excitation  coefficients.  The  waveguide  mode 
excitations  generally  are  correlated,  and  so  the  stochastic  average  value 
<  ?(  0) ,  6,  iji)>  of  the  power  density  is  written  in  the  form 


'■P(ij,e,y)>  =  I  ,>;  <a  ,  vw)  >  <a  (u)>P  , 
y  p  p  p  UP 


+  Z  a^((jj)P  (u,0,$) 

u  u  uu 


+  Z’^Z'c  ,  (<u)P  ,  (^,8.<j,) 
y'u  p'u 


(5-23) 


where 


(  (D  )  =  the  variance  of  the  y mode  coefficient, 

U 

C  foj)  “  the  covariance  function  for  the  mode  coefficient 

^  ^  and  the  conjugate  of  the  y  '  mode,  and 


where  the  primes  on  the  last  double  sum  signify  that  only  terms  y  ^  y' 
are  included  in  the  sum. 

Equations  (5-21)  and  (5-23)  are  the  appropriate  equations  for 
obtaining  average  power  pacterns  over  designated  frequency  intervals  for 
deterministic  and  rand\.?  «ituatioQ8,  respectively.  This  process  has 
been  described  and  numerically  simulated  in  Section  II  Cor  a  phased 
array  antenna.  Equations  for  time  domain  analysis  are  needed. 
Accordingly,  attention  is  focused  on  the  time  domain  response  in  the 
following  paragraphs. 

The  time  domain  electric  field  E(t,  6  ,^  )  can  be  computed  as  the 
discrete  Fourier  transform  of  the  frequency  domain  response 


E(  ,v 


,  $  )  as 


E(t,9,<»  =  —  L  E((i)  ,9,^)H(u  )exp[ju  t] 
iS  jj  n  n  n 


(5-24) 
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where  H(  “  )  is  the  pulse  spectrum  function  and  where  the  summation  is 
over  a  set  of  frequencies  exceeding  twice  the  bandwidth  of  the  input  pulse 
spectrum.  The  f(  9  ,  <!'  )  in  Equation  (5-24)  are  the  frequency  domain 
fields  obtained  when  the  input  frequency  spectrum  is  uniformly  weighted 
in  <wq>litude.  The  time  domain  power  density  P(t,  9.  'J')  proportional  to 
E*(t,  6  ,  ♦Itit,  e  ,  41)  and  is  written  explicitly  as 


P(t.6.4)  =-^2  E*{u.^..e.$>E(«^.6,4)H*(..^.)KCw^) 


expliCu)  -u)  ,)tj  }  (5-25) 

n  n 

This  equation  for  P(t,  0.  4)  can  be  expressed  directly  in  terms  of  the 
mode  propagation  coefficients  a(  u )  as 


P(t.e,*)  =  ^2  h  (0.4)  Z,l, 


(5-26; 


+  exp[jui  A“?(u.^..0..fr>A®(<o^.e.$) 

+  exp[-j»i)  ^”1  ^ 

(5-27) 
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and  where  all  other  symbols  have  been  previously  defined.  Equation  (S-26) 
is  the  appropriate  equation  for  computing  the  power  density  envelope  of 
the  radiated  pulse  at  the  far-field  observation  point  (  6|<|>  )  for  a 
specified  input  time  pulse  whose  spectrum  is  H(  ‘Ibe  pulse  amplitude 
envelope  is,  of  course,  the  square  root  of  P(t,6  ,  0).  The  statistical 
average  value  of  the  power  density  envelope  of  the  time  pulse  is  obtained 
by  Statistically  averaging  over  the  mode  coefficient  products  ,(ui^,) 
■  a (  111  ).  Making  use  of  the  fact  chat  the  mode  coefficients  for 
different  frequencies  are  uncorrelated,  the  statistical  average  value 
^  P(t,  6  ,iji  )>  of  Pit, 6  ,  4i)  is  written  as 


12  ,  * 

'P(t,0,4i)>  =  ;; — 2  h  (6,$)  Z  ,1  Z  ,1;  { <a  ,  (ui  ,)>  <a  (u  )> 
HTT  n  n  u  u  u  n  P  n 

(n’/  a) 


yt  1 

U)  ,U3  H  (uj  ,)H(a»  )exp{j(w  ,-u)  )r/c]*  P  ,  ,  (u)  ,  ,u  6,4i)  } 

n  n  n  n  n'  n  n  nw  p  n  n' 


+  Z  Z  ,Z  )H((u  ) 

n  u'  u  n  n  n 


-  p  ,  (ui.,w  c  ,  } 

n  DM  y  tf  n'  u  p 


(5-28) 


The  statistical  average  value  of  the  envelope  of  the  pulse  amplitude 
cannot  be  obtained  by  taking  the  square  root  of  Equation  (5-28)  because 
the  average  value  of  the  square  root  of  a  random  variable  is  not  equal 
to  the  square  root  of  the  average  value.  The  statistical  average  value 
of  the  amplitude  of  the  pulse  can  be  computed  numerically  with  the  aid 
of  the  Nakagami  p.d.f.  if  it  is  needed.  However,  the  statistical 
average  value  of  the  envelope  of  the  pulse  power  density  given  by 
Equation  (5-28)  is  the  quantity  of  interest  in  this  study. 

C.  Tism  Domain  Simuistions 

Numerical  computeciona  were  performed  to  compute  the  pulse 
envelope  pover  density  of  non-random  pulses  'sdiated  by  a  waveguide 
hom-fed  reflectoc  antenna  for  both  in-band  and  out-of-band  operation 
in  the  presence  of  a  nearby  metallic  cylinder.  The  comiiutations  vtte 
made  with  aid  of  existing  computer  programs  APATT  and  MBD  previously 
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developed  at  Georgia  Tech  I  6  ,38l.  APATT  was  used  to  generate  the 
dear-site  in-b<jnd  and  out-of-band  antenna  patterns  as  a  function  of 
frequency.  HBD  was  employed  to  compute  the  total  electric  field 
(antenna  ♦  scatterer)  versus  frequency  for  the  antenna/obstacle 

situations.  The  radiated  pulse  power  density  was  calculated  via 
Equation  (5-21). 

The  antenna  used  for  the  simulations  is  a  4-foot  diameter 

paraboloidal  reflector  antenna  fed  by  a  waveguide  horn.  The  design  in- 
band  center  rrequency  is  3.0  GHz  and  the  design  in-band  antenna 

polarization  is  vertical.  The  waveguide  portion  of  the  feed  system 

consists  of  a  6-inch  long  section  of  standard  WR-  S-band  waveguide 

excited  hy  a  coax-to-waveguide  adapter.  I'he  antenna  operates  in  the 
presence  of  a  12-foot  tall,  2-foot  diameter  metal  cylinder  located 

within  the  near-field  of  the  antenna. 

Plots  of  the  radiated  pulse  power  density  envelope  are  shown  in 
Figures  5-2  through  5-7  for  a  single  input  rectangular  pulse  of  two 
nanoseconds  duration.  Figures  5-2  and  5-3  show  the  radiated  pulses  for 
in-band  operation  of  the  antenna  while  Figures  !-4  through  5-7  show  the 

radiated  pulses  for  out-of-band  operation  of  the  antenna.  The  plots 

are  sequenced  as  foll<'ws.  Figures  5-2  and  5-3  show  Che  radiated  pulse 
in  the  antenna  boresight  (O-degree)  direction  for  the  clear  site  and 
Che  blocked  condition,  respectively.  Figures  5-4  and  5-5  show  the  out- 
of-band  radiated  pulse  in  the  boresight  direction  antenna  boresight 
for  dear-site  and  blocked  operation,  respectively.  Similarly,  Figures 
5-6  and  5-7  show  the  ouc-of-band  radiated  pulse  in  a  sidelobe  region 
for  the  clear  site  and  blocked  situations,  respectively.  For  the 
blocked  case,  the  metal  cylinder  was  located  six  feet  away  from  Che 
antenna  alcng  the  5°  radial.  The  four  out-of-band  figures  are 
sequenced  in  Che  same  manner  and  the  obstacle  was  located  at  the  sasm 
position  .iu3t  described  for  the  in-band  coa^mtations . 

The  in-band  pulses  were  computed  from  Che  antenna  or 
autenna/obstacle  electric  field  for  frequencies  from  2.0  GBx  to  4.0  GHz 
in  0.05  GHz  steps.  Similarly,  the  ouc-of-band  pulses  were  coaq>uCed 
from  the  antenna  or  antenna/obstacle  elecric  fields  for  frequencies 
froE  5.5cHzto7.5  GHr  in  0.05  GHz  steps.  The  pulse  carrier  frequencies 
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boresighi:  direction  for  the  unobsLi acted  4.0-foot  diameter 
reflector  entenna  for  an  input  2.u-nanoaecond  rectangular 
pulse  with  carrier  frequency  at  the  out-of-band  frequency 


CQP)d3M0d  3AIXV13U 
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fcr  the  in-band  and  out-of-band  pulses  are  3.0  CHs  and  6.5  CBa, 
respectively. 

The  waveguide  power  flow  for  Che  in-band  coaputaCions  was  entirely 
in  the  TE^q  node.  The  waveguide  power  flow  was  raoJonly  diatribuCed 
among  the  allowed  modes  for  the  out-of-band  computations.  In 
particular,  the  mode  coefficients  at  each  increawntal  frequency  were 
obtained  from  a  random  draw  via  a  pseudo-random  number  generator.  The 
statistical  average  power  for  each  allowed  mode  was  chosen  to  be 
(I.O/N)  watt,  where  N  is  the  total  cumber  of  allowed  modes  at  a 
particular  frequency. 

The  ill-band  pulses  displayed  in  Figures  5-2  and  5-3  are  noticeably 
distorted.  In  particular,  the  radiated  pulses  for  both  the  dear-site 
and  partially  blacked  situations  are  characterised  by  the  appearance  of 
time  side lobes  and  noticeable  increase  is  the  pulsewidth,  as  expected 
from  theory.  The  pulse  for  the  partially  blocked  antenna  suffers  less 
distortion  than  the  pulse  for  Che  unobstructed  antenna.  This  result  is 
attributed  to  the  fact  chat  Che  variations  in  the  joiplitudea  and  phases 
of  frequency  domair.  electric  fields  are  less  dramatic  for  the  partially 
blocked  antenna  than  for  the  unobstructed  antenna. 

The  out-of-band  pulses  displayed  in  Figures  5-4  through  5-7  are 
severely  distorted.  The  out-of-band  pulses  in  the  antenna  boreaight 
direction  are  considerably  more  distorted  than  their  in-band 
counterparts.  Of  course,  this  result  is  expected  because  the  electric 
field  varies  more  dramatically  over  the  out-of-band  frequency  interval 
than  for  Che  in-band  frequency  interval.  The  out-of-band  pulae  in  the 
bores ight  direction  for  Che  unobstructed  antenna  does  not  resemble  Che 
input  2.0  nanosecond  rectangular  pulse.  However,  the  out-of-band  pulse 
in  the  boresighC  direction  for  the  partially  blocked  antenna  does 
resemble  a  rectangular  pulse  and  Che  distortion  is  manifested  priiurily 
in  the  appearance  of  time  sidelobea.  The  tread  toward  smaller  pulse 
distortion  for  the  partially  blocked  antenna  for  the  pulsos  in  the 
boresight  direction  is  also  evident  in  Figures  5-6  and  5-7  for  Che 
pulses  radiated  in  the  sidelobe  directions. 

The  observed  trend  Coward  leas  distortion  for  Che  partially 
blocked  antenna  Chan  for  Che  unobstructed  sntenna  ie  attributed  Co  the 


moderating'  effect  of  the  scattering  obstacle  on  the  variations  of  the 
amplitude  and  phases  of  the  electric  field  over  the  2.0  GHz  frequency 
interval.  However,  it  should  be  noted  that  the  2.0  foot  cylinder  is 
large  ii\  terms  of  wavelength  and,  consequently,  the  scattered  fields 
are  less  frequency  sensitive  chan  the  scattered  fields  for  resonant 
scacterers.  Accordingly,  It  would  be  premature  to  generalise  the 
distortion  reaults  to  othei  scattering  environments  involving 
structures  that  have  diisenaiona  cotsparable  to  the  wavelengths  in  Che 
pulse  spectrum,  particularly  those  environments  where  dielectric 
scatterers  are  present.  Additional  analytical,  numerical,  and 
experimental  research  is  needed  in  order  to  characterize  the  wideband 
cw  and  pulsed  effects  of  arbitrary  scattering  environstents. 
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SECTIOK  VI 

CONCLUDING  REMARKS  AND  RECOMMENDATIONS 


The  research  work  described  herein  meets  or  exceeds  the  stated 
objectives  for  this  basic  research  program.  The  overall  suin  objective 
was  to  advance  the  state  of  knowledge  concerning  aear~£ield  theory  and 
techniques  for  wideband  radiating  r;ystems  at  in-band  and  ouc-of-bar.d 
frequencies.  This  was  achieved  by  the  successful  completion  of  four 
research  tasks.  The  task  objectives  and  the  key  results  for  each  task 
are  described  below.  These  results  are  followed  by  recesmsendat ions  for 
farther  research  work  in  the  four  task  areas. 

A.  Suasaary  of  Results 

Task  1.  Provide  a  near-field  methodology  to  characterize 

electromagnetic  emitter  radiation  patterns  at  in-band  and  out-of-band 
frequencies  for  wide  bandwidth  radiators.  The  objective  of  this  task 
was  to  develop  the  appropriate  theory  and  equations  based  on 
statistical  analysis  techniques  for  efficient  characterization  of 
wideband  radiators. 

The  theory  and  equations  were  developed  for  characterizing  the 
radiation  patterns  of  wideband  cw  or  pulsed  antennas  over  both  in-oand 
and  out-of-band  frequency  intervale  from  measured  data  collected  ”ia 
near-field  measurement  techniques.  The  results  are  applicable  to 
either  phased  array  or  reflector  antennas.  Numerical  simulations  were 
performed  for  (1)  a  20-element  out-of-band  waveguide  phased  array  with 
no  inter-elesient  coupling.  The  effect  of  statistical  correlations  of 
the  near-field  data  was  studied  and  methods  for  handling  correlation 
effects  were  derived.  Useful  approximations  for  the  probabi  lity 
density  function  for  the  radiated  power  pattern  statistics  for 
correlated  sources  were  also  identified. 

The  key  results  obtained  via  the  Task  1  efforts  may  be  stzmsarized 
as  follows: 

(i)  Tha  statistical  average  patterns  and  standard  deviations 
at  selected  frequencies  can  provide  a  very  succinct 
engineering  descrirtion  cf  the  important  SIC  characteristics 
of  wideband  CW  multimoding  antennas.  The  statistical 
average  patterns  and  standard  deviations  are  a  practical 


alCernaclve  to  the  comparatively  more  expensive  and 
cwsbersome  Monte  Carlo  simulations. 

(2)  the  statistical  average  pattern  for  a  given  frequency 
may  be  computed  from  a  knowledge  of  the  following  near~ 
field  statistical  parax^ters: 

(a)  statistical  average  value  of  the  electric  field 
at  all  near-field  measurement  points, 

(b)  the  standard  deviation  of  the  electric  field  at 
all  measurement  points,  and 

(c)  the  covariance  functions  for  the  electric  fields 
at  all  different  near-field  measuratent  points. 

(3)  The  ststistical  average  pattern  versus  time  for  a  pulsed 
system  depends  on  all  of  the  above  near-field  stastistical 
parameters  listed  ia  Item  2,  and  the  following  far-field 
statistical  parameters : 

(a)  the  statistical  average  value  of  the  far-field 
electric  field  at  all  frequencies  in  the  frequency 

bsnd ,  , 

(b)  the  standard  deviation  of  the  far-field  electric 
field  at  all  frequencies  in  the  frequency  band,  and 

(c)  the  covariance  functions  of  the  electric  field  at 
all  different  frequencies  in  the  frequency  band. 

However,  the  far-field  statistical  parameters  listed  above 
can  be  computed  from  the  near-field  statistical  paraiaeters. 

Thus,  pulsed  antennas  characterisation  does  not  require 
knowledge  of  any  ad:Zitional  statistical  near-field  data. 

(4)  The  probability  density  function  (p.d.f.)  for  the  far-field 
electric  fields  of  cor~elated  random  sources  or  scatterers 
can  be  accurately  approximsted  for  most  antenna  or  scattering 
problems  via  the  use  of  Makagami’s  p.d.I.  snd  linear  operator 
theory.  A  less  accurate  but  relatively  simple  expression  for 
the  p.d.f.  was  also  derived  that  is  suitable  for  certain 
engineering  applications. 

Task  2.  Theoretically  relate  the  radiation  pattern 
characterisation  to  the  basic  data  needed  for  efficient  optimisation 
analysis  of  electromagnetic  spectrum  usage.  Tbe  objective  of  this  task 
was  to  relate  the  near-field  derived  wideband  antenna  characterisation  » 

to  antenna  coupling  of  antenna  systeis  which  co-exist  in  the  same  BM 
environment .  , 
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The  efforts  on  this  task  %rere  devoted  to  deriving  approximate  but 
accurate  coupling  prediction  equations..  Three  different  techniques  for 
coupling  analysis  were  considered.  These  three  techniques  are  denoted 
respectively  as  (1)  Plane  Wave  Spectrum  (FWS),  (2)  Spherical  Spectrum 
Wave  (SWS) ,  and  (3)  the  Geometrical  Theory  of  Diffraction  (GTD).  Most 
of  Che  research  efforts  were  devoted  toward  development  of  the  PWS 
technique.  Multiple  scattering  effects  are  not  addressed  in  the 
analysis.  Bowever,  multiple  scattering  effects  can  be  approximately 
analysed  if  the  scattering  matrix  of  each  antenna  is  specified  from 
theory  or  measurements.  Theory  and  equations  for  the  FWS,  SWS,  and  GTD 
approaches  were  presented  and  discussed.  Results  of  numerical 
simulstions  using  the  PWS  approach  for  analyzing  antenna  coupling  are 
also  presented.  The  GTD  technique  for  deterministic  antenna  analysis 
was  adapted  to  yield  comparatively  sisiple  equations  for  certain  classes 
of  antenna  coupling  problems. 

The  analyses  for  this  tank  were  conducted  in  the  frequency  domain. 
The  priinary  goal  was  to  obtain  equations  valid  over  wide  frequency 
bandwidths .  The  tiaw  domain  response  can  then  be  obtained  by 
numerically  computing  the  Fourier  Transform  of  the  frequency  domain 
equations.  The  wideband  frequency  response  of  swept  Git  radiating 
aystems  is  then  characterized  by  numerically  computing  the  frequency- 
averaged  pattern  and  standard  deviation. 

Task  3.  Provide  the  methodology  to  assess  the  effects  of  system 
devices  (i.e.,  higher-order  mode  generation)  on  the  radiation  pattern. 
The  two  main  objectives  fur  this  task  trere  (1)  to  develop  a  method  for 
analyzing  the  antenna  pattern  effects  of  higher  order  modes  generated 
and/or  propagated  by  cossKn  feed  system  devices  snd  (2)  to  assess  the 
feasibility  of  deducing  the  mode  excitations  for  a  radiating  feed 
system  from  analysis  of  measured  out-of-band  pattern  data. 

Equations  were  derived  during  Task  1  and  Task  2  to  permit  the  out- 
of-band  radiation  pattern  statistics  to  be  computed  frc»  a  knowledge  of 
the  system  aevice  statistics  which  describe  the  higher-order  mode 
generation  and  propagation.  Theory  and  equations  for  identifying  the 
higher-order  mode  excitations  of  a  feed  containing  passive  and/or 
active  waveguide  devices  through  utilization  of  wideband  or  pulsed 
near-field  measurements  were  also  derived. 
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The  ciaely  coaipletion  of  the  required  research  work  for  this  task 
permitted  a  brief  digression  into  another  iaportant  aspect  of  the  out- 
of-band  response  of  systea  devices.  In  particular,  equations  were 
sought  to  describe  the  influence  of  particular  waveguide  devices  on  the 
excitation  and  subsequent  propagation  of  out-of-band  waveguide  aodes. 
Three  waveguide  devices  were  considered  in  this  exploratory  study, 
naotely  (1)  coax-to—waveguide  adapters,  (2)  radia)  bends,  and  (3)  ferrite 
phase  shifters.  Prograw  constraints  precluded  a  detailed  quantitative 
analysis  of  all  three  kinds  of  devices.  Accordingly,  the  research 
efforts  were  concentrated  mainly  on  the  theoretical  and  numerical 
analysis  of  a  coax-to—waveguide  adapter. 

Task  4.  Investigate  the  impact  of  site  effects  on  the  near-field 
antenna  analysis  technology,  The  objective  of  this  task  was  to  extend 
the  existing  moncchrosut ic  spectrv»  scattering  matrix  analysis  to  study 
antenna  siting  effects  on  the  wideband  and  out-of-band  performance  of 
radiating  systems. 

A  theoretical  study  was  perforved  to  investigate  techniques  and 
derive  the  equations  for  extending  the  existing  TOoochromatic  spectrum 
scattering  matrix  theory  to  analyse  site  effects  over  wide  frequency 
bandwidths.  In  particular,  frequency  doswin  equations  were  derived  via 
the  Plane  Wave  Spectrum  Scattering  (PVSS)  technique.  The  pulse 
envelope  response  is  obtained  froc  the  Fourier  Transform  of  the 
frequency  domain  WSS  equations.  Based  upon  the  results  of  the 
theoretical  study,  a  numerical  simulation  of  the  site  effects  on  the 
wideband  in-band  and  out-of-band  performance  of  a  pulsed  radiating 
system  was  performed  for  a  reflector  antenna  partially  blocked  by  a 
round  atetal  cylinder  located  in  the  antenna's  near-field. 

B.  Recosmendations 

The  results  of  the  research  work  conducted  for  this  basic  research 
program  significantly  advance  the  state  of  knowledge  concerning  near¬ 
field  wideband  in-band  and  out-of-band  radiation  and/or  acattering 
phenomena.  Howver,  further  advances  can  and  should  be  sude  through 
additional  research  work.  Accordingly,  recowtedatioas  for  further 
research  work  are  presented  herewith; 
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Task  1 


(a)  Conduct:  theoretical  and  numerical  analyses  to  derive  a 
reliable,  accurat<^  practical  method  to  compute  the  near¬ 
field  covariance  functions  from  the  near-field  data 
obtained  from  a  single  probe. 

(b)  Conduct  theoretical  and  numerical  analyses  to  derive  the 
probability  density  function  for  correlated  sources  or 
scatterers . 

(c)  Conduct  theoretical  and  numerical  analyses  to  derive  the 
inter-element  coupling  and  covariance  functions  for  a 
multimoding  array  of  out-of-band  waveguide  elements. 

(d)  Conduct  near-field  measurements  for  a  multimoding  waveguide 
array  to  provide  validation  data  for  item  (c). 

Task  2. 

(a)  Conduct  theoretical  and  numerical  analyses  to  derive  antenna 
coupling  model (s)  that  are  valid  for  arbitrary  orientations 
and  locations  of  in-band  as  out-of-band  wideband  cw  or  pulsed 
reflector  antennas  in  the  presense  of  environmental 
scattering  obstacles. 

(b)  Conduct  near-field  antenna  coupling  measurements  in  the 
presence  of  selected  scattering  obstacles  in  order  to 
provide  validation  data  for  item  (a). 

Task  3. 

(a)  Conduct  theoretical  and  numerical  analyses  to  device  improved 
approximate  methods  for  computing  the  higher-order  mode 
excitations  and  propagation  constants  for  common  waveguide 
components  such  as  radial  bends,  ferrite  phase  shifters, 
coax-to-waveguide  adapters,  and  waveguide  rotary  joints. 

(b)  Conduct  near-field  measurements  for  multimoding  feed  systems 
to  generate  validation  data  for  item  (a)  and  to  provide 
insight  into  realistic  mode  excitations  encountered  in 
real-world  feed  systems. 

Task  4. 

(a)  Conduct  theoretical  and  numerical  analyses  to  derive 
simplified  expressions  for  making  rapid  but  accurate 
approximate  analyses  of  the  wideband  cw  or  pulsed  response 
of  reflectors  or  phased  array  antennas  that  are  partially 
obstructed  by  one  or  more  near-field  obstacles. 
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(b)  Conduct  experisents  to  aeature  the  wideband  cw  or 
pulsed  response  of  partially  blocked  antennas  to 
provide  validation  data  for  ites  (a). 
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